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warping  was  established.  A  revelation  from  these  closed-form  solutions  is 
that,  elastic  coupling  lowers  tlie  first  coupled  frequencies  (in  fact  a 
significant  amount  of  coupling  could  reduce  the  first  frequency  to  almost 
zero.  The  discovery  of  the  closed-form  solutions  for  the  free  vibration 
which  seems  to  mark  the  first  time  such  solutions  were  ever  obtained,  not 
only  led  to  answers  to  a  number  of  previously  unanswered  questions  biit  also 
raised  new  unanswered  questions  sucli  as  "Does  the  aeroelastic  divergence  or 
flutter  problem  of  such  wings  have  any  closed-form  solutions". 

During  tite  second  year  of  the  first  phase  investigation  was  concerned  with 
essentially  how  to  find  the  answer  to  the  above  mentioned  question.  Ac¬ 
cordingly  in  the  first  quarter, 
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efforts  were  concentrated  on  determining  the  possibility  of  finding  closed-form  solutions  to 
the  divergence  problem.  The  investigation  led  to  two  possible  methods  of  obtaining  such 
solutions.  These  methods  are  (a)  the  "elimination"  approach  (which  was  used  for  the  free 
vibration)  and  (b)  the  method  of  Laplace  transformation.  Although  it's  already  known  that  (a) 
works,  (b)  was  implemented  in  formulating  the  analytical  expressions  for  the  closed-form 
solution  of  the  divergence  problem  due  to  an  anticipated  relative  ease.  In  the  beginning  of  the 
2nd  phase,  the  closed  form  divergence  and  vibrations  results  which  have  and  are  still  being 
extracted  are  being  studied  very  carefully  to  establish  physical  trends  in  the  constraint  of 
warping  phenomena.  In  particular,  the  computed  results  are  being  studied  in  order  to  evolve 
the  basic  mechanism  by  which  the  constrained  warping  model  may  be  used  to  explain  the 
concept  of  aeroelastic  tailoring.  Preliminary  results  that  are  being  evolved  seem  to  indicate 
that  in  free  vibration,  if  the  elastic  coupling  is  varied,  there  are  some  circumstances  when  the 
first  two  coupled  modes  may  diverge  from  one  another  and  other  circumstances  when  they 
may  approach  coalescence.  The  big  question  therefore  is  "would  such  coalescence  result  in 
instability"  even  in  absence  of  air  (or  at  low  air  speeds)?  Another  question  whose  answer  is 
being  searched  also  is  "what  is  the  role  of  such  coalescence  or  diverging  potential  in 
aeroelastic  tailoring"?  The  answers  to  these  questions  seem  to  be  evolving  gradually  as  the 
investigation  continues.  For  example,  it  is  seen  as  shown  in  this  report  that  coalescence  of 
modes  may  not  result  in  stabilities  (even  at  very  low  flight  velocities). 

in  the  most  recent  investigations,  efforts  have  been  concentrated  on  the  unsteady 
aerodynamics.  This  was  considered  necessary  especially  since  this  research  program  plans 
to  eventually  look  at  the  aeroelastic  behavior  of  wings  in  transonic  flow  -  a  non-linear,  very 
complicated  flow  regime.  It  is  recognized  that  problematic  phenomena  at  low  subsonic 
speeds  can  potentially  be  disastrous  at  transonic  speeds.  Hence  a  new  discovery  by  the 


principal  investigator  in  non-linear  fluid  dynannics  is  employed  to  establish  some  closed  form 
solutions  for  2-D  unsteady  non-linear  flow  applicable  to  transonic  regime  which  are  expected 
to  be  eventually  used  for  aeroelastic  analysis.  Some  of  the  results  have  been  compiled  and 
published  in  the  AIAA  Journal.  In  the  process  of  studying  the  unsteady  2-D  nonlinear 
transonic  flows,  what  seems  to  be  a  major  breakthrough  was  discovered:  this  breakthrough 
is  the  discovery  that  the  3-D  steady  non-linear  transonic  flow  equations  can  be  transformed 
into  a  linear  hodograph  equivalent.  This  discovery  essentially  reverses  a  hundred  years  state 
of  thinking  and  belief  in  the  scientific  and  mathematical  community  that  the  hodograph 
method  was  limited  to  the  2-D  flows.  The  results  of  this  discovery  have  been  written  up  and 
published  in  the  AIAA  Journal. 

In  addition  another  paper  AIAA  92-257  which  presents  practical  shock  free  3-D  wing  for 
transonic  flight  is  to  be  presented  in  June,  1992  at  the  AIAA  applied  Aerodynamics 
Conference  in  Palo  Alto,  California.  This  research  program  seems  to  be  seeing  the  beginning 
of  what  seems  to  be  another  breakthrough,  perhaps  even  more  significant  than  what  we  have 
seen  thus  far.  This  is  what  seem  to  be  a  discovery  that  transformations  methods  can  even 
be  used  to  effectively  study  the  ultimate  set  of  equations  in  continuum  gas  dynamics  known 
as  the  Navier  Stokes  set  of  equations.  This  new  finding  could  be  the  key  to  effectively 
unlocking  the  secret  of  gas  dynamics  and  aeroelastic  phenomena  for  real  gases  with  viscosity, 
which  have  been  blurry  at  best  thus  far.  This  finding  has  resulted  in  a  preliminary  paper 
accepted  to  be  presented  at  an  international  conference  to  be  held  in  Colorado  in  August  of 
1992.  It  is  being  proposed  that  the  investigation  of  these  new  findings  be  pursued  in  a 
P''oposed  next  phase  of  research.  This  is  because  if  the  findings  are  correct  this  could  be  the 
start  of  a  new  era  of  fluid  or  gas  dynamics  analysis. 


1 .0  ABSTRACT 


The  investigation  of  an  aeroelastically  induced  constrained  warping  phenomenon  for 
a  composite,  supermaneuverable  type  aircraft  wing  has  continued  in  this  second  phase  of  the 
study.  The  first  phase  investigation  was  concentrated  mainly  on  the  static  phenomena  and 
the  search  for  closed  form  solutions  for  free  vibration  of  aircraft  wings  having  constrained 
warping  in  the  presence  of  elastic  coupling.  The  wing  is  analytically  modelled  as  a  straight 
flat  laminated  plate.  Various  forms  of  highly  simplified  aerodynamic  loads  are  employed  in  the 
analysis.  The  free  vibrations  and  stability  aspects  of  this  phenomenon  are  examined  to  obtain 
some  physical  insights  and  to  determine  its  importance  and/or  design  implications.  Analytical 
tools  employed  include  an  affine  transformation  concept  which  was  formulated  previously  (by 
the  present  principal  investigator)  as  well  as  a  non-dimensionalization  scheme.  With  the  help 
of  these  tools,  an  evolution  of  effective  warping  parameters  with  which  to  study  this 
phenomenon  was  carried  out.  The  virtual  work  theorem  and  variational  principles  were  used 
to  derive  the  equations  of  motion  based  on  the  assumed  wing  displacements.  The  existence 
of  closed-form  free  vibrations  solutions  for  composite  wings  with  elastic  coupling  and 
constraint  of  warping  was  established.  A  revelation  from  these  closed-form  solutions  is  that, 
elastic  coupling  lowers  the  first  coupled  frequencies  (in  fact  a  significant  amount  of  coupling 
could  reduce  the  first  frequency  to  almost  zero.  The  discovery  of  the  closed-form  solutions 
for  the  free  vibration  which  seems  to  mark  the  first  time  such  solutions  were  every  obtained, 
not  only  led  to  answers  to  a  number  of  previously  unanswered  questions  but  also  raised  new 
unanswered  questions  such  as  "Does  the  aeroelastic  divergence  or  flutter  problem  of  such 
wings  have  any  closed-form  solutions." 

During  the  second  year  of  the  first  phase,  investigation  was  concerned  with  essentially 
how  to  find  the  answer  to  the  above  mentioned  question.  Accordingly  in  the  first  quarter. 


efforts  were  concentrated  on  determining  the  possibility  of  finding  closed-form  solutions  to 
the  divergence  problem.  The  investigation  led  to  two  possible  methods  of  obtaining  such 
solutions.  These  methods  are  (a)  the  "elimination"  approach  (which  was  used  for  the  free 
vibration)  and  (b)  the  method  of  Laplace  transformation.  Although  it's  already  known  that  (a) 
works,  (b)  was  implemented  in  formulating  the  analytical  expressions  for  the  closed-form 
solution  of  the  divergence  problem  due  to  an  anticipated  relative  ease.  In  the  beginning  of  the 
2nd  phase,  the  closed  form  divergence  and  vibrations  results  which  have  and  are  still  being 
extracted  are  being  studied  very  carefully  to  establish  physical  trends  in  the  constraint  of 
warping  phenomena.  In  particular,  the  computed  results  are  being  studied  in  order  to  evolve 
the  basic  mechanism  by  which  the  constrained  warping  model  may  be  used  to  explain  the 
concept  of  aeroelastic  tailoring.  Preliminary  results  that  are  being  evolved  seem  to  indicate 
that  in  free  vibration,  if  the  elastic  coupling  is  varied,  there  are  some  circumstances  when  the 
first  two  coupled  modes  may  diverge  from  one  another  and  other  circumstances  when  they 
may  approach  coalescence.  The  big  question  therefore  is  "would  such  coalescence  result  in 
instability"  even  in  absence  of  air  (or  at  low  air  speeds)?  Another  question  whose  answer  is 
being  searched  also  is  "what  is  the  role  of  such  coalescence  or  diverging  potential  in 
aeroelastic  tailoring"?  The  answers  to  these  questions  seem  to  be  evolving  gradually  as  the 
investigation  continues.  For  example,  it  is  seen  as  shown  in  this  report  that  coalescence  of 
modes  may  not  result  in  stabilities  (even  at  very  low  flight  velocities). 

In  the  most  recent  investigations,  efforts  have  been  concentrated  on  the  unsteady 
aerodynamics.  This  was  considered  necessary  especially  since  this  research  program  plans 
to  eventually  look  at  the  aeroelastic  behavior  of  wings  in  transonic  flow  -  a  non-linear,  very 
complicated  flow  regime.  It  is  recognized  that  problematic  phenomena  at  low  subsonic 
speeds  can  potentially  be  disastrous  at  transonic  speeds.  Hence  a  new  discovery  by  the 
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principal  investigator  in  non-linear  fluid  dynamics  is  employed  to  establish  some  closed  form 
solutions  for  2-D  unsteady  non-linear  flow  applicable  to  transonic  regime  which  are  expected 
to  be  eventually  used  for  aeroelastic  analysis.  Some  of  the  results  have  been  compiled  and 
Dublished  in  the  AIAA  Journal.  In  the  process  of  studying  the  unsteady  2-D  nonlinear 
transonic  flows,  what  seems  to  be  a  major  breakthrough  was  discovered:  this  breakthrough 
is  the  discovery  that  the  3-D  steady  non-linear  transonic  flow  equations  can  be  transformed 
into  a  linear  hodograph  equivalent.  This  discovery  essentially  reverses  a  hundred  years  state 
of  thinking  and  belief  in  the  scientific  and  mathematical  community  that  the  hodograph 
method  was  limited  to  the  2-D  flows.  The  results  of  this  discovery  have  been  written  up  and 
published  in  the  AIAA  Journal. 

In  addition  another  paper  AIAA  92-257  which  presents  practical  shock  free  3-D  wing  for 
transonic  flight  is  to  be  presented  in  June,  1992  at  the  AIAA  applied  Aerodynamics 
Conference  in  Palo  Alto,  California.  This  research  program  seems  to  be  seeing  the  beginning 
of  what  seems  to  be  another  breakthrough,  perhaps  even  more  significant  than  what  we  have 
seen  thus  far.  This  is  what  seem  to  be  a  discovery  that  transformations  methods  can  even 
be  used  to  effectively  study  the  ultimate  set  of  equations  in  continuum  gas  dynamics  known 
as  the  Navier  Stokes  set  of  equations.  This  new  finding  could  be  the  key  to  effectively 
unlocking  the  secret  of  gas  dynamics  and  aeroelastic  phenomena  for  real  gases  with  viscosity, 
which  have  been  blurry  at  best  thus  far.  This  finding  has  resulted  in  a  preliminary  paper 
accepted  to  be  presented  at  an  international  conference  to  be  held  in  Colorado  in  August  of 
1992.  It  is  being  proposed  that  the  investigation  of  these  new  findings  be  pursued  in  a 
proposed  next  phase  of  research.  This  is  because  if  the  findings  are  correct  this  could  be  the 
start  of  a  new  era  of  fluid  or  gas  dynamics  analysis. 
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=  chordwise  integrals 

=  affine  space  half-chord  and  chord,  respectively 
=  elastic  constants 

=  parameter  that  measures  the  location  of  the 
reference  axis  relative  to  mid -chord 

=  bending  and  torsional  stiffness,  respectively 
=  wing  box  depth 

=  affine  space  bending  and  torsional  displacement, 
respectively 

=  elastic  coupling  and  warping  stiffness,  respectively 
=  elemental  stiffness  parameter 
=  Strouhal  number 

=  affine  space  running  aerodynamic  lift  and  moments, 
respectively 

=  affine  space  half-span  for  the  wing 
=  affine  space  mass  per  unit  span 

=  differential  aerodynamic  pressure  distributions  in 
physical  and  affine  space,  respectively 


t 


=  time 


U.Uo  =  virtual  work  expressions  in  physical  and  affine  space, 

respectively 

U,  =  flutter  speed 

(x,y,z),  (Xo,yo,Zo  =  physical  and  affine  space  coordinates,  respectively 

Y  ,,  P ,  =  displacement  shape  functions 

r,L„L2,D*,D*o,  =  generic  nondimensionalized  stiffness  parameters 
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=  affine  mass  ratio  parameter 

=  Poisson  ratios  and  generalized  Poisson's  ratio, 
respectively 

=  divergence  parameter 

=  affine  space  material  and  air  density,  respectively 
=  twisting  displacements 
=  displacement 
=  vibration  frequency 


3.0  INTRODUCTION 


In  the  late  1 960's  designers  began  to  investigate  the  possibility  of  exploiting  the  directional 
properties  of  composite  materials  to  improve  aeroelastic  characteristics  of  lifting  surfaces.  For 
example,  superfighter  designers  may  aeroelastically  tailor  a  wing  so  that  it  deforms  to  the 
optimum  camber  under  maneuvering  loads. 

For  instance,  a  wing  can  be  tailed  so  that  its  leading  edge  will  twist  downward  under  the 
stress  of  a  tight  turn,  thereby  decreasing  the  wing's  angle  of  attack  and  hence  reduce  drag. 

Exploitation  of  the  directional  properties  of  composite  materials  to  solve  aeroelastic 
instability  problems  (known  as  aeroelastic  tailoring)  may  enhance  the  performance  of  future 
high  performance  military  aircraft,  particularly  the  super-maneuverable  concepts.  This  stems 
from  the  need  for  future  military  weapon  systems  to  exhibit  high  performance  and  minimum 
vulnerability  (e.g.:  minimum  radar  cross-sectional  area).  These  requirements  may  lead  to 
unorthodox  aircraft  configurations  which  in  turn  result  in  unorthodox  aeroelastic  instability 
problems  (e.g.,  the  X-29  primary  mode  of  instability,  called  the  'body  freedom'  flutter). 

Perhaps  the  aeroelastic  tailing  concept  would  have  been  discovered  much  earlier  if  the 
physics  of  anisotropic  aeroelasticity  were  more  apparent.  This  inherent  physical  intractability 
is  largely  due  to  the  existence  of  numerous  variables,  e.g.,:  flight  parameters,  several 
composite  directional  properties,  fiber  orientation  angles,  etc.,  which  are  not  even  truly 
independent,  in  the  anisotropic  system.  This  state  of  affairs  is  therefore  analogous  to  that 
which  existed  in  the  basic  rigid  body  aerodynamics  before  the  advent  of  the  similarity  rule 
theory.  This  theory  clearly  revealed  that,  the  utilization  of  non-dimensionalized  groups  of 
variables,  e.g.:  Reynolds,  Mach,  Strouhal,  Froude  numbers  etc.,  provide  significantly  superior 
physical  insight  to  the  problem  than  the  individual  physical  variables.  The  new  methodology 
that  is  being  used  as  the  basic  tool  in  this  research  program  is  basically  the  aeroelastic 


equivalent  of  the  aerodynamic  similarity  rule.  The  expected  superiority  of  this  new  approach 
over  the  state-of-the-art  (SOA)  counterpart,  which  utilizes  individual  physical  variables, 
especially  in  terms  of  physical  insights,  has  been  demonstrated  in  references  1-10.  For 
example,  a  high-aspect-ratio  composite  wing  could  behave  aeroelastically  like  a  low  aspect 
ratio  wing  and  vice-versa.  Similarity  parameters  can  expose  conditions  for  which  this  might 
happen.  This  is  significant,  (for  instance)  in  the  light  of  the  important  role  played  by  the  wing 
aspect  ratio  in  the  aerodynamics  approximations  for  an  aeroelastic  analysis.  This  result  may 
therefore  be  suggesting  some  new  form  of  coupling  between  the  elastic  and  aerodynamic 
equations  in  composite  wing  aeroelasticity.  A  fundamental  aspect  of  this  new  methodology 
has  been  used  in  studies  at  Purdue  University  ”  and  MIT’^. 

In  this  research,  an  investigation  of  a  wing's  spanwise  sectional  distortion  (warping) 
resulting  from  aerodynamic  forces  and  its  influence  on  the  wing's  free  vibrations,  as  well  as 
(aeroelastic)  stability,  is  being  carried  out.  In  particular,  St.  Venant's  torsional/twisting  theory, 
which  is  currently  widely  used  for  estimating  the  wing's  twisting  displacements  has  been 
examined  with  the  help  of  the  new  methodology,  to  determine  its  limitations,  when  applied 
to  wings  fabricated  of  composite  materials.  The  relevance  and  significance  of  this  study  for 
the  newly  emerging  supermaneuverable  type  aircraft  may  be  seen  in  the  light  of  the  fact  that 
(a)  supermaneuverability  is  characterized  by  high  angle  of  attack  which  implies  high  twisting 
aerodynamic  forces,  and  (b)  most  aircraft  designers  believe  future  aircraft  will  be  fabricated 
of  40-70%  composite  materials. 

When  the  St.  Venant's  torsional  theory  is  used  to  estimate  a  wing's  twisting  displacement 
and/or  forces,  the  fact  that  the  wing's  root  section's  distortion  is  relatively  small  compared 
to  that  of  other  sections  is  ignored.  However,  previous  investigators  have  determined  that 
such  an  unrealistic  assumption  may  lead  to  only  little  errors  if  the  aspect  ratio  of  the  wing  is 
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very  high. 

The  research  has  shown  during  the  first  phase  that  the  conclusions  reached  by  previous 
investigators  are  basically  true  for  wings  fabricated  of  metals  or  isotropic  materials.  A  set  of 
new  theories  is  therefore  being  postulated  in  this  research  effort  for  accurately  estimating  the 
twisting  displacements,  vibrational  frequencies  and  instability  boundaries  for  wings  fabricated 
of  composite  materials. 

In  the  second  phase,  the  mechanism  of  instability  is  defined  as  the  main  target  of  more 
detailed  investigation.  This  is  partially  because  some  new  findings  during  this  research 
program  seem  to  indicate  that  composite  wings  undergoing  free  aeroelastic  oscillations  exhibit 
some  unique  phenomena  not  known  to  occur  in  metal  wings.  One  of  this  phenomena  is  what 
we  have  referred  to  as  a  "damping  like  behavior"  (in  absence  of  actual  damping)  brought 
about  by  the  presence  of  asymmetric  coupling.  This  among  other  things  presents  opportunity 
for  a  merging  of  two  modes  during  a  free  vibration.  The  natural  question  is  therefore  "would 
such  a  phenomenon  lead  to  instabilities"  as  has  been  well  documented  for  other  aeroelastic 
systems?  Of  primary  interest  to  an  aeroelastician  is  "how  does  this  phenomenon  connected 
to  aeroelastic  tailoring?  These  are  some  of  the  questions  that  we  are  searching  for  answers 
to  in  this  second  phase  of  this  research.  It  is  beginning  to  be  apparently  evident  that  modal 
coalescence  is  possible.  Furthermore  it  is  also  evident  that  this  modal  coalescence  may  lead 
to  some  form  of  instability  at  low  flight  velocities  which  we  are  still  trying  to  discern  properly 
and  carefully. 

In  preparation  for  more  sophisticated  flutter  or  divergence  analysis,  a  new  discovery  by  the 
principal  investigator  in  nonlinear  unsteady  fluid  dynamics  is  being  refined  to  enhance  and  help 
expose  the  physical  nature  of  aeroelastic  instability  of  composite  wings  with  elastic  coupling 
and  warping  restraint.  Some  of  the  preliminary  results  have  just  been  published  in  the  AIAA 


Journal.  A  copy  of  this  publication  is  included  in  this  writeup.  Furthermore,  the  study  of  the 
unsteady  transonic  fluid  dynamics  has  led  us  to  some  new  findings.  In  the  process  of 
studying  the  unsteady  2-D  nonlinear  transonic  flows,  what  seems  to  be  a  major  breakthrough 
was  discovered:  this  breakthrough  is  the  discovery  that  the  3-D  steady  non-linear  transonic 
flow  equations  can  be  transformed  into  a  linear  hodograph  equivalent.  This  discovery 
essentially  reverses  a  hundred  years  state  of  thinking  and  belief  in  the  scientific  and 
mathematical  community  that  the  hodograph  method  was  limited  to  the  2-D  flows.  The 
results  of  this  discovery  have  been  written  up  and  published  in  the  AIAA  Journal. 

In  addition  another  paper  AIAA  92-257  which  presents  practical  shock  free  3-D  wing  for 
transonic  flight  is  to  be  presented  in  June,  1992  at  the  AIAA  applied  Aerodynamics 
Conference  in  Palo  Alto,  California.  This  research  program  seems  to  be  seeing  the  beginning 
of  what  seems  to  be  another  breakthrough,  perhaps  even  more  significant  than  what  we  have 
seen  thus  far.  This  is  what  seem  to  be  a  discovery  that  transformations  methods  can  even 
be  used  to  effectively  study  the  ultimate  set  of  equations  in  continuum  gas  dynamics  known 
as  the  Navier  Stokes  set  of  equations.  This  new  finding  could  be  the  key  to  effectively 
unlocking  the  secret  of  gas  dynamics  and  aeroelastic  phenomena  for  real  gases  with  viscosity, 
which  have  been  blurry  at  best  thus  far.  This  finding  has  resulted  in  a  preliminary  paper 
accepted  to  be  presented  at  an  international  conference  to  be  held  in  Colorado  in  August  of 
1992.  It  is  being  proposed  that  the  investigation  of  these  new  findings  be  pursued  in  a 
proposed  next  phase  of  research.  This  is  because  if  the  findings  are  correct  this  could  be  the 
start  of  a  new  era  of  fluid  or  gas  dynamics  analysis. 
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II 


4.0  RESEARCH  OBJECTIVES 


The  overall  goal  of  this  phase  of  the  research  program  is  twofold:  first,  is  to  formulate 
closed'form  solutions  to  the  aeroelastic  divergence  and  flutter  of  aircraft  wings  with  warping 
constraint  and  elastic  coupling,  and  use  these  to  study  the  effects  of  constrained  warping  on 
aeroelastic  response;  second  is  to  develop  new  plate-beam  finite  elements  and  use  in  modal 
analysis  of  vibration  and  flutter  to  provide  an  independent  check  on  the  closed-form  solution 
results. 

During  a  preliminary  investigation  of  the  constrained  warping  under  aeroelastic  forces, 
a  phenomenon  being  studied  in  the  current  research  program  sponsored  by  the  AFOSR,  the 
Principal  Investigator  discovered  (for  the  first  time),  that  there  are  closed-form  solutions  to  the 
free  vibrational  and  divergence  problem  of  composite  aircraft  wings  having  warping 
constraints  and  elastic  coupling  [2J.  In  this  study  these  newly  discovered  techniques  are  to 
be  extended  (and  expanded)  and  used  for  continuing  the  formulation  of  closed-form  solutions 
to  the  aeroelastic  flutter  and  dynamic  response  problems  of  aircraft  wings  with  the  constraint 
of  warping  and  elastic  coupling.  This  process  is  expected  to  involve  determining  the 
necessary  symmetries  that  exist  in  the  flutter  and  dynamic  response  equations  (by  comparing 
them  to  their  free  vibrations  and  divergence  counterpart)  to  permit  the  extraction  of  their 
characteristic  roots  for  elastically  coupled  deformations  of  the  wing  under  the  influence  of 
aerodynamic  forces. 

Experience  from  the  free  vibrational  and  the  ongoing  divergence  investigations  has 
shown  that  the  computations  involved  in  extracting  the  eigenvalues  from  the  closed-form 
eigen  expressions  derived  using  the  new  techniques  are  very  challenging.  For  example,  it  is 
difficult  to  determine  how  to  eliminate  or  minimize  the  numerical  errors  resulting  from  the 
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prcfo/^re  -f'^rihe  -f  ^n/ 

of  dynamic  response  problems  are  to  be  evolved  with  the  help  of  the  experience  gained  from  the 
study  of  the  vibrations  and  divergence.  The  closed-form  solutions  thus  extracted,  are  to  be 
thoroughly  studied  to  understand  the  physics  of  the  constrained  warping  phenomena.  In 
particular,  the  computed  results  used  to  determine  the  basic  mechanism  by  which  the  constrained 
warping  model  can  be  used  to  explain  the  concept  of  aeroelastic  tailoring.  This  is  necessary 
because  the  basic  mechanism  of  aeroelastic  tailoring  seems  to  be  still  rather  elusive.  For 
example,  it  seems  to  be  generally  accepted  that  tailoring  occurs  because  ply  angles  are  "properly” 
oriented.  From  our  ongoing  research,  tailoring  seems  to  be  more  than  just  a  ply  angle 
orientation.  Careful  study  of  the  data  generated  are  expected  to  be  carried  out  to  try  to  unfold 
this  apparent  mystery.  In  other  words,  by  using  a  figure  similar  to  Fig.  (1),  how  can  tailoring 
(higher  flutter  speed)  be  explained? 


<^10.  <^20 '  isotropic  1st  & 
2nd  frequency,  respectfully. 

<^11.  <«>22  -  composite  1st  & 
2nd  frequency,  respectfully. 

Vpi  =  flutter  speed 


13 


5.0  STATUS  OF  RESEARCH  EFFORTS 
5.1  INTRODUCTION 

During  the  reporting  period,  the  research  program  progressed  according  to  plan  resulting  in 
the  accomplishment  of  the  goals  defined  for  the  period.  These  goals,  which  have  been  defined 
in  Section  4  of  this  report,  basically  are  to  extend  and  utilize  some  of  the  newly  developed 
important  tools  in  the  current  ongoing  research  program  which  is  being  sponsored  by  AFOSR 
to  investigate  the  constrained  warping  phenomena  of  aircraft  wings  fabricated  of  composite 
materials  having  elastic  coupling  and  subjected  to  aeroelastic  forces.  The  proposed  study  is 
basically  comprised  of  the  investigation  of  (a)  the  basic  mechanism  of  instability  and  dynamic 
response,  (b)  the  mechanism  of  modal  transformation  in  flutter  and  dynamic  response,  and 
(c)  nonlinear  and  transonic  instabilities. 

In  the  first  and  second  reporting  period  during  the  preliminary  study.  Prof.  G.  A.  Oyibo 
started  the  investigation  of  the  basic  mechanism  of  aeroelastic  instability  of  composite  aircraft 
wings  with  restrained  warping  and  elastic  coupling.  Following  a  realization  that  proper 
understanding  of  the  mechanism  of  aeroelastic  instability  particularly  the  dynamic  instabilities, 
require  some  sound  analytical  tools,  a  review  of  the  important  tools  developed  for  the 
divergence  and  free  vibration  for  composite  wings  with  warping  restraint  and  elastic  coupling 
during  the  last  phase  was  started.  During  this  review,  a  thorough  study  of  the  preliminary 
data  generated  were  carried  out.  The  understanding  gained  during  this  review  provided  the 
necessary  background  necessary  for  properly  formulating  the  instability  mechanism  problem. 
Consequently  an  expansion  or  extension  of  these  earlier  tools  started.  First,  the  major 
scientific  revelation  to  the  scientific  community  made  possible  (through  the  help  of  these 
tools)  by  this  research  program  was  properly  reexamined.  The  revelation  is  that  a  composite 
wing  having  a  restrained  warping 


in  the  presence  of  an  asymmetric  elastic  coupling  (in  absence  of  damping)  exhibit  a  "damping¬ 
like"  behavior.  This  new  discovery  which  should  be  critical  to  the  proper  understanding  of  the 
mechanism  of  aeroelastic  tailoring  is  still  being  thoroughly  investigated.  The  latest  results 
show  that  modal  coalescence  is  possible.  These  results  are  seen  to  show  that  such  a 
coalescence  could  lead  to  instability  (even  at  low  flight  speeds).  More  careful  investigations 
are  now  being  carried  out  in  order  to  "fine  tune”  the  conclusions  of  these  results  and  to 
recommend  possible  experimental  studies  to  verify  them.  In  addition,  a  refinement  of  a  new 
discovery  in  nonlinear/transonic  fluid  dynamics  by  the  principal  investigator  was  initiated  in 
anticipation  of  more  realistic  or  sophisticated  fiutter/divergence  analyses  to  further  the  goals 
of  this  research.  Some  of  the  preliminary  results  have  been  compiled  and  published  in  the 
AIAA  Journal.  In  the  process  of  studying  the  unsteady  2-0  nonlinear  transonic  flows,  what 
seems  to  be  a  major  breakthrough  was  discovered:  this  breakthrough  is  the  discovery  that 
the  3-D  steady  non-linear  transonic  flow  equations  can  be  transformed  into  a  linear  hodograph 
equivalent.  This  discovery  essentially  reverses  a  hundred  years  state  of  thinking  and  belief 
in  the  scientific  and  mathematical  community  that  the  hodograph  method  was  limited  to  the 
2-D  flows.  The  results  of  this  discovery  have  been  written  up  and  published  in  the  AIAA 
Journal.  In  addition  another  paper  AIAA  92-257  which  presents  practical  shock  free  3-D  wing 
for  transonic  flight  is  to  be  presented  in  June,  1992  at  the  AIAA  applied  Aerodynamics 
Conference  in  Palo  Alto,  California.  This  research  program  seems  to  be  seeing  the  beginning 
of  what  seems  to  be  another  breakthrough,  perhaps  even  more  significant  than  what  we  have 
seen  thus  far.  This  is  what  seem  to  be  a  discovery  that  transformations  methods  can  even 
be  used  to  effectively  study  the  ultimate  set  of  equations  in  continuum  gas  dynamics  known 
as  the  Navier  Stokes  set  of  equations.  This  new  finding  could  be  the  key  to  effectively 
unlocking  the  secret  of  gas  dynamics  and  aeroelastic  phenomena  for  real  gases  with  viscosity. 
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which  have  been  blurry  at  best  thus  far.  This  finding  has  resulted  in  a  preliminary  paper 
accepted  to  be  presented  at  an  international  conference  to  be  held  in  Colorado  in  August  of 
1992.  It  is  being  proposed  that  the  investigation  of  these  new  findings  be  pursued  in  a 
proposed  next  phase  of  research.  This  is  because  if  the  findings  are  correct  this  could  be  the 
start  of  a  new  era  of  fluid  or  gas  dynamics  analysis. 

The  investigation  of  the  numerical  method  for  efficiently  and  accurately  extracting  the 
eigenvalues  for  the  vibration  and  divergence  problems  continued.  New  numerical  methods  for 
handling  some  problem  regions  of  these  problems  are  also  being  examined.  Prof.  T.  A. 
Weisshaar  is  also  continuing  his  investigations  into  how  a  new  beam-plate  finite  elements  may 
be  formulated  to  capture  the  constraint  of  warping  effect  in  order  to  provide  results  which 
may  be  used  to  independently  check  the  closed-form  solutions  which  have  been,  and  are 
being  generated. 


5.2  ACCURATE  DIVERGENCE  THEORY  FOR  COMPOSITE  SUPERMANEUVERABLE 
AIRCRAFT  WINGS 
5.2.1  INTRODUCTION 

Modern  supermaneuverable  aircraft  concepts  benefit  a  great  deal 
from,  among  other  things,  significant  advances  in  materials 
technology  and  the  availability  of  more  accurate  aerodynamic  pre¬ 
diction  capabilities.  Supermaneuverability  as  a  design  goal  in¬ 
variably  calls  for  an  optimization  of  the  design  parameters. 
Optimization  may  be  partially  accomplished  for  example,  by  using 
composite  materials  to  minimize  weight.  Indeed,  it  has  been  known 
that  these  composite  materials  can  be  tailored  properly  to  resolve 
the  dynamic  or  static  instability  problems  of  these  types  of  air¬ 
craft.  The  concept  is  referred  to  as  aeroelastic  tailoring. 

While  aeroelastic  tailoring  has  tremendous  advantages  in  the  design 
of  an  aircraft,  the  analysis  which  provides  the  basis  for  the 
aeroelastic  tailoring  itself  is  generally  very  involved.  This,  is 
rather  unfortunate  since  a  good  fundamental  physical  insight  of 
the  tailoring  mechanism  is  required  for  accurate  and  reliable  re¬ 
sults  . 

In  this  investigation  an  attempt  is  made  to  look  at  some  dyneunics 
theories  that  can  be  used  to  understand  the  aeroelastic  tailoring 
mechanism.  Specifically,  the  accuracy  of  the  St.  Venant  torsion 
theory  which  is  relatively  simple  and  frequently  used  in  aeroelas¬ 
tic  analysis  is  examined  with  particular  reference  to  the  effects 
of  the  wings  aspect  ratio  as  well  as  other  design  parameters. 
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An  accurate  torsion/twist  theory  is  particularly  significant  for 
supermaneuverable  aircraft  wings  since  supermaneuverability  is 
basically  characterized  by  high  angle  of  attack. 

Although  earlier  studies  (1,2,3)  have  indicated  that  the  St.  Ve- 
nant's  torsion  theory  is  reasonably  accurate  except  for  aircraft 
wings  with  fairly  low  aspect  ratios,  the  theory  supporting  that 
conclusion  was  based  on  the  assumption  that  the  wing  is  construc¬ 
ted  of  isotropic  materials.  Basically  the  St.  Venant's  torsion 
theory  assumes  that  the  rate  of  change  of  the  wing's  twist  angle 
with  respect  to  the  spanwise  axis  is  constant.  This  assumption 
is  hardly  accurate  particularly  for  modern  aircraft  construction 
in  which  different  construction  materials  are  employed  and  the 
aerodynamic  loads  vary  significantly  along  the  wing's  span.  How¬ 
ever,  References  1-3  have  shown  that  (in  spite  of  such  an  inac¬ 
curate  assumption)  the  main  parameter  that  determines  the  accur¬ 
acy  of  the  St.  Venant's  theory  is  the  wing's  aspect  ratio.  Thus, 
it  was  determined  that  the  theory  is  fairly  accurate  for  moderate 
to  high  aspect  ratio  wings  constructed  of  isotropic  materials. 

In  recent  studies  (4,5,6)  however,  it  has  been  shown  that  for 
wings  constructed  of  orthotropic  composite  materials,  the  con¬ 
clusion  of  References  1-3  need  to  be  modified.  Rather  than  using 
the  geometric  aspect  ratio  of  the  wing  to  determine  the  accuracy 
of  St.  Venant's  twist  theory,  it  was  suggested  that  a  generic 
stiffness  ratio  as  well  as  an  effective  aspect  ratio  which  consid¬ 
ers  the  wing's  geometry  and  the  ratio  of  the  principal  directional 
stiffness  should  be  considered  in  establishing  the  accuracy  of  St. 


Venant's  theory. 


The  present  investigation  is  related  to  the  studies  that  were  ini¬ 
tiated  in  References  5  and  6.  In  this  study  the  first  task  was  to 
examine  the  role  of  coupling  (both  mass  and  elastic  coupling)  on 
the  accuracy  of  St.  Venant's  theory  applied  to  static  problems. 

It  was  discovered  that  coupling  plays  a  very  significant  role  on 
the  accuracy  of  St.  Venant's  twist  theory.  The  second  task  was  to 
investigate  the  torsional  vibration  for  a  flat  plate  model  of  an 
aircraft  wing  fabricated  of  composite  materials  in  which  the  con¬ 
strained  warping  phenomenon  is  more  realistically  represented, 
with  particular  emphasis  on  higher  frequencies  and  to  compare  re¬ 
sults  with  those  from  a  representation  based  on  St.  Venant's 
theory . 

5.2.2  FORMULATION 

Consider  an  aircraft  wing  fabricated  of  composite  materials  and 
mathematically  idealized  as  a  cantilevered  plate  subjected  to  an 
aerodynamic  flow  over  its  surfaces.  The  mathematical  statement  of 
the  virtual  work  theorem  for  such  a  plate  model  is  well  known  and 
documented.  It  is  also  known  that  such  mathematical  statements 
of  the  virtual  work  theorem  for  a  laminated  plate  model  are  char¬ 
acterized  with  the  existence  of  so  many  variables  (in  the  state¬ 
ment)  ,  reflecting  the  various  directional  properties  for  the  leunin- 
ated  plate  model,  which  would  tend  to  interfere  with  any  physical 
insight  that  might  be  desired  from  a  phenomenological  analysis 
employing  such  a  mathematical  statement.  The  newly  discovered  af¬ 
fine  transformation  concept  (5,6  and  7)  was  developed  principally 
to  resolve  such  a  problem. 


If 


This  new  concept  therefore  can  be  used  to  evolve  the  mathematical 
statement  of  the  virtual  work  theorem  in  an  affine  space  given  by 
the  following  equation. 


are  the  elastic  constants,  p,  is  the  material  density,  Ap 
is  the  differential  pressure  distribution,  w  is  the  displacement, 
t  is  the  time,  A  integrals  represent  area  integrals  and  h  is  the 
wing  box  depth. 

Equations  (1)  and  (2)  therefore  form  the  basis  of  the  newly  de¬ 
veloped  methodology.  The  equations  of  motion  of  a  plate  model  of 
an  aircraft  wing  can  now  be  derived  by  prescribing  a  realistic 
wing  displacement  and  using  Equation  (1) . 

When  Equation  (1)  is  compared  to  its  physical  space  counterpart, 
it  is  seen  that  Equation  (1)  has  fewer  variables.  It  is  also  seen 
that  Equation  (1)  contains  only  non-dimensionalized  stiffness  quan¬ 
tities  (compared  to  dimensional  stiffness  quantities  in  its  physi¬ 
cal  space  counterpart) .  Another  feature  of  this  new  methodology 
which  makes  it  unique  is  that  the  non-dimensionalization  (a  conse¬ 
quence  of  the  affine  transformation)  is  accomplished  before  assum¬ 
ing  the  wing  deformations.  This  means  that  the  non-dimensionaliza¬ 
tion  is  independent  of  how  the  wing  deforms.  A  non-dimensionaliza¬ 
tion  scheme  that  depends  on  a  particular  assumption  of  the  wings 
deformations  could  lead  the  analyst  to  an  incorrect  physical  inter¬ 
pretation  of  results,  since  the  wing's  deformations  assumptions 
have  inherent  errors  because  they  are  based  on  the  analyst's  judg¬ 
ments  and  experience.  This  observation  may  become  clearer  during 
the  evolution  of  a  warping  parameter  with  which  to  study  the  ef¬ 
fects  of  the  warping  constraint  phenomenon  on  the  status  and  dyna¬ 
mics  of  a  wing  fabricated  of  composite  materials  later  in  this  study. 
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If  the  chordwise  curvature  is  neglected  in  an  initial  approxima¬ 
tion,  the  wing's  deformation  may  be  assumed  as  follows: 


w(t,x^,y^) 


X  a  (t,y  ) 
0  o'  ■'o' 


(3) 


where  h^  and  are  the  bending  and  twisting  displacements,  re¬ 
spectively. 


It  can  be  shown  that  when  Equation  (3)  is  substituted  into  Equation 
(1)  and  the  variational  calculus  is  carried  out  for  arbitrary 
and  a^,  the  following  equations  of  motion  are  obtained: 


a.h  a,  «'"+  P„a,h  +  p„a,a 

I  /%  7  rt  h  rt  I  rt  •  Ck  / 


1v 


iii, 


■l"o 


’2  0 


5“o  "  ^o“l"o  ■  ^o“2“o 


^3"^  ‘  U^o'  *  Po"3®o  ^  Po^z\  “  % 


(4) 


where : 


Z2 


5.2.3 


EVOLUTION  OF  WARPING  PARAMETERS 


The  evolution  of  the  warping  parameter  with  which  to  study  the  con¬ 
strained  aeroelastic  warping  phenomenon  for  wings  fabricated  of 
composite  materials  is  a  process  that  depends  on  the  sophistication 
of  the  wing's  mathematical  model;  whether  coupling  effects  are  in¬ 
cluded,  whether  the  wing's  chordwise  curvatures  are  included  and 
so  on.  Therefore,  any  warping  parameter  is  as  good  as  the  corres¬ 
ponding  wing's  displacements  assumptions.  However,  Equation  (1) 
makes  it  possible  for  the  analyst  to  determine  its  effective  inde¬ 
pendent  variables  even  before  the  displacement  assumptions  are  made. 


By  non-dimensionalizing  the  spanwise  space  variable  in  Equation  (4) , 
depending  on  whether  one  is  interested  in  the  static,  dynamic, 
coupled  or  uncoupled  displacements,  one  of  the  following  warping 
parameters  may  be  useful. 


■f^o 


(7) 


-A 

C  ■  Cg  vj 


3  /  n*  _  ^2  \ 


(8) 


A.  =  — 
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3  A. 

1-332 


(9) 


•v..' 
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where : 


0„  =  0  (1  .«) 


(10) 


(So/Cq)  i®  defined  as  the  wing's  effective  aspect  ratio  and  D^* 
and  L  are  the  generalized  stiffness  and  coupling  ratios,  respec¬ 
tively  (defined  in  earlier  work  such  as  References  5  and  6) . 


Equations  (7)  thru  (9)  represent  the  appropriate  warping  parame¬ 
ter  for  dyncunic  deformation,  static  displacement  with  elastic 
cross-coupling,  and  static  deformation  with  "geometric"  coupling 
{e  1)  ,  respectively. 


It  was  discovered  in  this  study  that  evolving  the  warping  parame 
ter  in  a  manner  shown  in  Equations  (1)  thru  (3) ,  should  enable 
one  to  investigate  the  effects  of  warping  on  the  composite  wing' 
dynamics  (or  the  accuracy  of  St.  Venant's  theory)  effectively. 
From  the  lamination  theory  for  composites  it  is  known  that  while 
D  *  and  (£  /c  )  are  always  positive,  L  and  a_  can  be  positive  or 
negative.  However,  from  Equations  (1)  and  (3),  it  is  clear  that 
whether  a  composite  wing  has  positive  or  negative  coupling,  the 
warping  effect  (in  terms  of  X^)  is  unchanged. 


Using  the  Laplace  transform  method  on  the  divergence  form  of  equa¬ 
tions  4  in  which  strip  aerodynamic  theory  is  used,  the  following 
Laplace  function  becomes  important  in  the  divergence  problem, 
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Equation  50  is  therefore  expected  to  rep  sent  the  closed-form  ex¬ 
pression  for  the  divergence  eigenvalues.  merical  methods  are  now 
being  used  to  extract  ti  ?se  eigenvalues  in  order  to  obtain  the 
divergence  speeds. 


RESULTS 


The  results  obtained  so  far  seem  to  verify  existing  results  and 
establish  new  trends.  Due  to  the  fact  that  the  evolution  of 
these  new  trends  are  yet  to  be  fully  completed,  only  a  summary 
of  these  results  are  presented  below. 

The  results  obtained  by  solving  equation  50  were  carefully  to 
evolve  the  physical  insight  into  the  mechanism  of  divergence  in¬ 
stability  in  the  presence  of  warping  restraint  and  elastic  coup¬ 
ling.  This  study  revealed  the  following,  as  can  be  seen  from 
figures  1-7, 

First,  it  is  seen  that  the  view  held  by  many  analysts  that  elas¬ 
tic  coupling  plays  a  significant  role  in  aeroelastic  divergence 
tailoring  is  verified. 

Second,  it  is  seen  that  another  view  that,  higher  aeroelastic 
divergence  stability  boundaries  are  feasible  with  negative  elas¬ 
tic  coupling  (than  positive  elastic  coupling) ,  is  basically  true, 
but  UP  to  a  point.  It  is  further  seen  that  there  seems  to  be  a 
limit  to  how  negative  the  elastic  coupling  can  be  made  to  obtain 
better  stability  boundaries  -  after  such  a  limit,  a  further  nega¬ 
tive  increment  of  elastic  coupling  would  seem  to  lower  the  stabil¬ 
ity  boundaries. 

Third,  it  is  found  that  the  effective  aspect  ratio  defined  in  the 
first  phase  of  this  research  program  (^^)»  for  simpler  models  can 
still  be  used  in  this  relatively  more  complex  model,  to  measure 
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the  effect  of  warping  restraint  on  the  phenomenon  of  divergence 
instability.  The  results  show  that  ignoring  the  warping  re¬ 
straint  would  lead  to  conservative  estimates  for  the  divergence 
instability  boundaries.  It  is  also  seen  that  the  restraint  of 
warping  effects  are  more  significant  for  small  effective  aspect 
ratio  (A^)  and/or  large  elastic  coupling. 
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Figure  4  Variation  of  divergence  dynamic  pressure  with  coupling 
for  large  effective  wing  aspect  ratio 
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Figure  5  Variation  of  divergence  dynamic  pressure  with  the  wings 
effective  aspect  ratio 
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Figure  6  Variation  of  divergence  dynamic  pressure  with  coupling 
for  large  effective  wing  aspect  ratio 
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Figure  7: 


Variation  of  the  nondimensionalized  dynamic 
pressure  with  wings  effective  aspect  ratio 
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Exact  Solutions  to  Aeroelastic  Oscillations  of  Composite  Aircraft 
Wings  with  Warping  Constraint  and  Elastic  Coupling^ 

h  * 

Gabriel  A.  Oyibo  and  James  Bentson 
Polytechnic  University 
Farmingdale,  New  York 

Abstract 

Exact  solutions  within  the  fraimework  of  standard  aeroelastic 
bending  and  twisting  assumptions  are  found  to  the  free  oscilla¬ 
tions  of  composite  aircraft  wings  having  warping  constraint  and 
elastic  coupling.  The  problem  is  treated  as  a  regular  boundary 
value  problem  consisting  of  two  fourth  order  partial  differen¬ 
tial  equations  coupled  by  the  presence  of  elastic  coupling. 

This  system,  which  is  linear,  therefore  is  equivalent  to  an 
eightii order  ordinary  differential  equation  .  Classical  linear 
"operator"  method  is  therefore  used  to  extract  fundeunental  solu¬ 
tions  which  are  superimposed  appropriately  to  obtain  an  exact 
functional  form  for  the  mode  shapes.  These  mode  shapes  are 
therefore  made  to  satisfy  the  necessary  boundary  conditions,  a 
process  that  leads  to  the  formulation  of  the  required  eigenvalue 
problem.  The  eigenvalues  are  extracted  numerically  by  using  ap¬ 
propriate  ordering  of  the  eight  roots  of  the  operator  equation. 
The  bending-torsion  frequencies  obtained  as  a  result  of  this  an¬ 
alysis  are  compared  favorably  with  existing  results.  New  in¬ 
sights  made  possible  by  these  results  which  are  preliminary,  ap- 


^Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research 
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pear  to  be  that  (a)  the  first  coupled  frequency  decreases  with 
increasing  coupling  and  (b)  the  phenomenon  of  modal  transforma¬ 
tions  found  by  earlier  investigators  is  explainable  in  terms  of 
some  conservative  inter-modal  energy  transfer. 


5.3.1  NOMENCLATURE 


a . 
1 


D.  . 
ID 


e 


=  chordwise  integrals 

=  affine  space  half-chord  and  chord,  respectively 
=  elastic  constants 

=  parameter  that  measures  the  location  of  the 
reference  axis  relative  to  mid-chord 


El,  GJ 
h 


K,So 

k.  . 
ID 


o  o 


m 


(Ap,Ap^) 


t 

U,U, 


=  bending  and  torsional  stiffness,  respectively 

=  wing  box  depth 

=  affine  space  bending  and  torsional  displacement, 
respectively 

=  elastic  coupling  and  warping  stiffness,  respec¬ 
tively 

=  elemental  stiffness  parauneter 

=  Strouhal  number 

=  affine  space  running  aerodynaunic  lift  and  mo¬ 
ments,  respectively 

=  affine  space  half-span  for  the  wing 

=  affine  space  mass  per  unit  span 

=  differential  aerodynamic  pressure  distributions 
in  physical  and  affine  space,  respectively 

=  time 

=  virtual  work  expressions  in  physical  and  affine 
space,  respectively 


U 


f 


(x,y,z,),  ^^o'^o^^o^ 


r  ,Lj^, L2  , D*  ,  D*^, 


u 


o 


X 

o 


flutter  speed 

physical  and  affine  space  coordinates,  respec¬ 
tively 

displacement  shape  functions 

generic  nondimensionalized  stiffness  pareuneters 
affine  mass  ratio  parameter 

Poisson  ratios  and  generalized  Poisson's  ratio, 
respectively 

divergence  parameter 


affine  space  material  and  air  density, 
spectively 

twisting  displacements 
displacement 
vibration  frequency 


5.3.2 

Introduction 


Perhaps  one  of  the  more  elusive  aspects  of  supermaneuverabil¬ 
ity  as  a  design  concept  is  its  aeroelastic  implications.  One  gen¬ 
erally  accepted  definition  of  supermaneuverable  aircraft  is  that  it 
is  designed  to  operate  at  high  angles  of  attack.  Strictly  speaking 
high  angle  of  attack  problems  are  nonlinear.  However  due  to  the 
high  degree  of  complexities  involved  in  dealing  with  nonlinear  aero¬ 
elastic  problems/ an  average  aeroelastician  would  prefer  to  deal  with 
a  linearized  version  of  the  problem  (at  least  as  a  first  approxima¬ 
tion)  .  If  linear  aeroelastic  equations  are  used  under  such  condi¬ 
tions,  at  least  it  should  be  assumed  that  the  high  angle  of  attack 
would  introduce  large  twisting  displacements  which  would  imply  that 
terms  containing  twisting  displacement  should  be  retained.  Even 
under  low  angle  of  attack  assumptions,  the  early  works  of  Reissner 
and  Stein^  and  later  works  of  Libescu  et  al^  have  shown  that  for 
metal  wings  there  are  conditions  under  which  the  so  called  St.  Ve- 
nant's  torsion  principle  is  inapplicable.  This  is  when  the  res¬ 
traint  of  warping  effect  is  important  and  a  more  accurate  analysis 
would  need  to  include  a  higher  order  term  involving  the  twisting 
displacement.  Although  the  retention  of  such  a  term  implies  solv¬ 
ing  fourth  order  (instead  of  second  order)  differential  equations 
for  the  twisting  and  bending  displacements,  the  equations  can  be 
easily  decoupled  for  metal  wings.  However  for  composite  wings,  the 
decoupling  of  these  equations  is  neither  easy  practically  nor  is  it 
even  desirable  from  aeroelastic  tailoring  standpoint.^  ^  These  stu¬ 
dies  have  also  shown  that  the  restraint  of  warping  is  very  important 


in  composite  wings.  Therefore  it  would  seem  that  an  aeroelastic 
analysis  of  a  supermaneuverable  {high  angle  of  attack)  aircraft 
wing  fabricated  of  composite  materials  would  need  to  consider 
the  effects  of  restraint  of  warping  as  well  as  elastic  coup¬ 
ling  . 

Previous  investigation  of  this  latter  problem  (free  vibra- 
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tion)  at  MIT  '  used  analytical  methods  to  solve  the  decoupled 
problem,  while  numerical  methods  were  utilized  to  solve  the  cou¬ 
pled  problems.  Consequently  general  results  were  presented  for 
the  decoupled  problem  while  representative  results  were  presen¬ 
ted  for  the  coupled  problem. 

In  this  paper  the  coupled  free  vibration  is  treated  analyti¬ 
cally  as  a  pair  of  coupled  fourth  order  (differential  equations) 
boundary  value  problem  to  which  exact  closed  form  eigen-solutions 
are  sought.  The  enforcement  of  the  necessary  boundary  conditions 
resulted  in  a  fairly  complicated  transcendental  function  to  be 
used  to  determine  the  required  eingenvalues  from  which  the  natur¬ 
al  frequencies  are  to  be  obtained.  This  transcendental  function 
was  complex  in  contrast  to  its  decoupled  counterpart  (which  is 
real) .  That  should  be  indicating  the  presence  of  the  phase  angle 
that  exists  between  the  twisting  and  bending  displacements.  A  com¬ 
parison  with  a  damped  (decoupled)  system  in  which  complex  determin¬ 
ant  signifies  phase  angles  between  damping  and  other  forces,  led  us 
to  the  formulation  of  an  explanation  for  the  "modal  transformation" 
phenomenon  which  was  reported  in  studies  at  MIT^''*  and  Purdue®  (which 
seemed  to  have  lacked  explanation  until  this  study) .  The  explana¬ 
tion  is  that  the  modal  transformation  may  be  viewed  as  a 
form  of  steady  state  conservative  (energy  stays  in  this  system 
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since  there  is  no  damping)  inter- 

modal  energy  transfer  between  the  vibration  modes.  In  fact 

0 

work,  currently  in  progress  at  Purdue  seem  to  support  and  con¬ 
firm  this  explanation.  The  results  which  favorably  compared 
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with  those  obtained  at  MIT,  '  also  revealed  that  coupling  has 
a  tendency  to  lower  the  first  coupled  natural  frequency  of  a 
composite  aircraft  wing.  In  fact  it  is  seen  that  a  substantial 
amount  of  coupling  could  reduce  the  first  coupled  frequency  to 
almost  zero  (hence  a  possibility  of  coupling  with  rigid  body 

modes)  . 

5.3.3 

Problem  Statement 

For  a  composite  aircraft  wing  cantilevered  at  the  root  as 
shown  in  Figure  1,  the  virtual  work  theorem  in  the  physical 
space  is  given  by 


6U  =  0  = 


§yyj'ki'‘"xx''  ^ 

0  A  ^ 


^°12  '''xx  "'yy 


°22^'''yy^  j 
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+  ^^16  '^'xy  ^^26  '''yy 

dxdydt  ~  I  J*  J  J  P  ^  w^dxdydt 


(1) 


•/// 

0  A 


0  A 


Ap6w  dxdydt 


where : 

D  are  the  elastic  constants,  p,  is  the  material  density,  Ap 

ij 

is  the  differential  pressure  distribution,  w  is  the  displacement, 
t  is  the  time,  A  integrals  represent  area  integrals  and  h  is  the 

wing  box  depth. 

Using  the  following  affine  transformation  of  variables; 


Dll  1/4 

X  = 

^22 


y  =  y. 


z  =  (2a) 


S3 


then  in  affine  space  the  virtual  work  theorem  becomes: 
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If  the  affine  space  equivalent  of  the  standard  aeroelastic 
displacement  assumptions  is  made,  i.e.. 


w(Xo,yo,t)  =  h^(y^,t)  +  x^a^(y^,t)  .. 


where  h  and  a  are  the  bending  and  twisting  displacements  re- 
o  o 

spectively,  then  it  can  be  shown  through  the  use  of  the  calcu¬ 
lus  of  variations  that  a  coupled  set  of  aeroelastic  equations 
of  motion  for  the  composite  aircraft  wing,  in  which  the  re¬ 
straint  of  warping  and  elastic  coupling  effects  are  accounted 
for  is  given  by 
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with  boundary  conditions 
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For  free  vibrations,  if  (through  the  geometric  construction 
of  the  wing)  is  made  to  be  zero,  equations  5  reduce  to 


i_  iv  ^  A  *  .  \  V 

lo  5o  olo  o 


(8a) 


iV  ,  ^  ^  L  **  U 
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with  boundary  conditions 


at  y  =  0, 
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a,h  =  ar-a  =  0,  aih 
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-  a^a'  =  0 


a,a^  +  aj-h^  -  a,a„  =  0 
i  O  3  0  4  o 


(8b) 


It  may  be  stated  here  that  the  restraint  of  warping  effect  is  rep- 

resented  by  the  product  of  a  and  a^  while  the  elastic  coupling 

effect  is  represented  by  the  a^  terms. 

5.3.4 

Methods  of  Solution 

Two  methods  for  solving  equations  8  are  examined  in  this  pa¬ 
per.  These  are  (a)  an  "exact  closed  form"  approach  and  (b)  a  "semi- 
exact  closed  form  approach.  The  exact  closed  form  approach  is  de¬ 
fined  here  as  one  in  which  explicit  expressions  are  derived  for 
the  eight  roots  of  the  eighth  order  operator  equations  represent¬ 
ing  equations  8a,  and  through  the  super imposition  of  fundamental 
solutions  corresponding  to  each  of  the  eight  operator  roots,  the 
boundary  conditions  8b  are  satisfied.  The  semi-exact  closed  form  ap¬ 
proach  is  the  same  as  (a)  except  that  the  roots  of  the  operator 
equation  are  determined  numerically  through  the  usage  of  some 
standard  root  extraction  subroutines. 

In  either  case,  to  solve  for  the  operator  roots  of  equation 
8a,  it  may  be  rewritten  in  operator  form  as  follows. 


and  Xjj  can  be  obtained  either  from  the  numerical  solution  of 
equation  11  or  according  to  the  method  described  in  Abramo- 


witz  and  Stegun,  as  follows 
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where  u.  (i  =  1,2,3)  is  the  root  that  makes  ^  r  c  all  real. 


5.3.5 

Consistency  Conditions 

From  equations  12(a)  and  12(b)  it  is  seen  that  there  are  six¬ 
teen  arbitrary  integration  constants  as  opposed  to  the  expected 
eight  constants.  The  additional  eight  constants  have  been  intro¬ 
duced  superf luausly  as  a  result  of  the  differential  operation 
which  was  done  in  order  to  eliminate  one  of  the  dependent  varia¬ 
bles  in  the  two  coupled  differential  equations  (equations  9) .  In 
order  to  get  rid  of  these  superfluous  solutions  it  is  necessary  to 
enforce  some  consistency  conditions.  This  may  be  accomplished  by 
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substituting  equations  12(a)  and  12 (b)  into  one  of  the  equations 
9  and  requiring  the  equation  to  be  satisfied  identically.  This 
procedure  shall  establish  a  set  of  explicit  relationships  be¬ 
tween  the  constants  and  in  which  can  be  determined  in 

terms  of  B  or  vice-versa, 
n 

When  the  superfluous  solutions  are  eliminated,  equations  12 
can  now  be  used  to  satisfy  the  boundary  conditions  for  this  prob¬ 
lem  (equations  8b) .  Consequently  the  condition  for  nontrivial  so¬ 
lutions  is  enforced  to  obtain  the  transcendental  functional  ex¬ 
pression  for  determining  the  eigenvalues  of  this  problem. 

5.3.6 

Eigenvalues 

The  following  steps  and  definitions  are  carried  out  in  order 
to  obtain  the  transcendental  functional  expression  for  determining 
the  necessary  eigenvalue  (for  this  eighth  order  boundary  value 
problem)  from  which  the  coupled  natural  frequencies  may  be  ob¬ 
tained 
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F53  =  -tj^a^^cosh3j^-t2a35Cos62+t^coshB3 
^63  =  tia26sinhP^-t2a4gSin62+t3sinhB3 
F^3  =  t3^aj^2CoshB3^+t2a32COsB2-t^cosB4 

^83  =  V28®^‘’^®1"4^48®^"V’^4®^"®4 
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Then,  are  given  by  the  roots  of 


F  = 


^^53^57'*'^63^67‘*'^73^  ^^54®58'’'^64®68'‘‘^84^ " 
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Equation  26  is  therefore  the  exact  closed  form  transcendental 
functional  expression  from  which  the  eigenvalues  may  be  ex¬ 
tracted.  The  coupled  natural  frequencies  are  related  to  these 
eigenvalues  by  the  following  expression. 


D„  1/2  u 
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5.3.7 

Computations 

The  extraction  of  the  eigenvalues  from  the  exact-closed  form 
transcendental  expression  in  equation  26  proved  to  be  a  very  chal 
lenging  computational  exercise  due  mainly  to  its  complex  nature, 
the  existence  of  branch  points,  the  necessity  to  order  the  opera 
tor  roots  appropriately  and  the  existence  of  numerical  noise. 
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The  characteristic  roots  for  the  equations  of  motion,  x^, 
were  extracted  in  two  different  ways  in  order  to  assure  accuracy: 
one  method  was  by  using  a  Jenkins-Traub  computational  method  from 
the  IMSL  library.  The  second  method  was  by  using  an  exact-closed 
form  approach  outlined  in  Abramowitz  and  Stegun.^^  When  these 
roots  are  being  extracted  numerically,  they  do  not  necessarily 
come  out  in  a  continuous  manner.  Therefore  a  subroutine  that  re¬ 
orders  them  so  that  they  become  continuous  with  respect  to  the 
eigenvalue  k^  is  also  employed.  Once  this  reordering  of  these 
roots  is  completed,  the  parameters  needed  for  computing  the  trans¬ 
cendental  function,  F,  in  equation  26  are  computed.  Finally,  the 
transcendental  expression  itself  is  computed  and  the  roots,  k^, 
are  found  numerically. 

One  interesting  result  is  that  the  values  of  F  are  complex 
here.  It  is  known,  however,  that  each  of  the  two  uncoupled  prob¬ 
lems  (bending  or  torsion)  when  treated  separately  has  a  real 
transcendental  expression  for  extracting  the  eigenvalues.  This  be¬ 
havior  of  the  transcendental  expression  for  extracting  the  coupled 
eigenvalues  (i.e.,  being  complex  as  opposed  to  being  real)  was  the 
first  hint  that  led  the  first  author  to  examine  any  possible  math¬ 
ematical  similarity  between  the  problem  at  hand  (a  coupled  non- 
damped  oscillations  problem)  and  simple  damped  oscillations  prob¬ 
lem  in  which  the  expression  for  determining  the  eigenvalues  is  in 
general  complex  due  to  the  need  to  determine  the  oscillation  fre¬ 
quency  and  the  amount  of  damping  in  the  mode,  etc.  In  the  coupled 
problem  the  complex  nature  of  this  transcendental  expression  for 
extracting  the  eigenvalues  seems  to  be  basically  a  reflection  of 


the  phase  that  normally  should  exist  between  the  bending  and  tor¬ 
sional  modes.  Realizing  that  in  the  damped  case,  there  is  a  non¬ 
conservative  energy  transfer  between  the  oscillating  system  and 
its  environment  as  well  as  intermodal  energy  transfer,  it  became 
clear  to  the  first  author  that  in  the  undamped  coupled  system 
there  may  be  a  conservative  energy  transfer  in  the  oscillating 
system.  At  the  time  of  this  thought  it  seemed  to  the  first  auth¬ 
or  that  if  it  made  sense,  the  idea  should  provide  an  explanation 

for  the  modal  changes  which  were  noticed  during  the  studies  at 
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Purdue  and  MIT  '  as  coupling  in  the  system  was  changed.  This 
prediction  which  turns  out  to  be  useful  and  confirmed  by  other 

g 

recent  studies  from  Purdue  University  shall  be  discussed  in  more 
detail  below. 

The  most  convenient  way  to  obtain  the  eigenvalues,  k^,  was 
found  to  be  by  graphical  means.  Thus  the  values  of  the  real  and 
imaginary  parts  of  the  complex  transcendental  function  F  are 
plotted  against  the  values  of  on  the  same  curve  as  shown  in 
Figure  2.  The  values  of  at  which  the  real  and  imaginary  parts 
of  the  transcendental  function  are  simultaneously  equal  to  zero 
corresponds  to  the  desired  eigenvalues  for  this  coupled  problem. 

One  of  the  problems  with  the  computational  model  described 
above  is  that,  when  the  coupling  parameter,  becomes  identical¬ 

ly  zero,  the  coupled  system  of  equations  becomes  computationally 
ill  conditioned  and  unsolvable.  To  circumvent  this  the  results 
for  the  uncoupled  case  ~  obtained  from  a  series  of 

calculations  using  successively  smaller  values  for  the  coupling 
parameter,  This  led  us  to  accept  the  values  of  eigenvalues 


for  zero  coupling  as  the  value  corresponding  to  the  limit  as  the 

coupling  approaches  zero.  Although  this  continuity  assumption 

seems  to  make  s'  -'se,  it  is  not  backed  by  a  rigorous  mathematical 

proof.  Luckily  it  was  possible  to  check  these  results  with  re- 

3  4 

suits  generated  for  an  isotropic/metal  aluminum/wing  by  MIT,  ' 
and  the  agreement  was  found  to  be  good  for  the  cases  checked. 

5.3.8 

Results  and  Discussions 

The  natural  frequencies  oi^  for  the  coupled  bending- tors  ion 

oscillations  for  a  composite  aircraft  wing  in  the  presence  of 

elastic  coupling  and  warping  restraint  is  found  (as  shown  by 

equation  27)  to  be  a  function  of  the  ratio  '  the  length 

i  and  the  nondimensionalized  frequency  parameter  k  for  the 
o  “ 

wing.  In  this  problem,  k^^,  is  a  function  of  only  two  parameters, 

i.e.,  X  which  may  be  considered  as  an  effective  nondimensional- 
o 

ized  aspect  ratio  and  which  in  a  way,  measures  the  amount 

of  elastic  coupling  in  the  wing  (X^/X^  =  1  for  zero  coupling) . 

It  is  therefore  seen  from  equation  27  that  in  order  to  increase 

w  one  needs  to  make  1  as  small  as  possible  and/or  make  (D22/P) 
n  ° 

and  k  as  large  as  possible.  Such  an  exercise  may  be  necessary 
n 

when  a  tailoring  of  the  frequency  is  needed  to  avoid  instabili¬ 
ties  (e.g.,  very  low  structural  frequencies  may  provide  an  atmos¬ 
phere  for  a  coupling  between  the  flexible  modes  and  rigid  body 
motions,  which  in  turn  has  a  potential  to  result  in  instability). 
This  kind  of  tailoring  is  made  convenient  through  the  use  of  equa¬ 
tion  27  in  which,  for  a  particular  wing  configuration  and  compos¬ 
ite  material,  every  variable  in  the  equation  shall  be  known  ex- 
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as  we  said  earli- 


cept  and  ^22’  Obviously  if  we  want  high  oj^, 

er,  0^2  should  be  made  to  be  as  high  as  possible  (and  of  course  p 

as  low  as  possible) .  Once  this  is  done  the  only  other  parameter 

to  be  tailored  is  k  . 

n 

The  plot  of  kj^  as  a  function  of  and  is  shown  in  Figure 

3  for  all  configurations  having  low  camber  to  twisting  coupling, 

and  all  values  of  bending  to  twisting  coupling.  The  results  com- 
3  4 

puted  at  MIT,  '  which  were  used  to  verify  the  present  results, 
are  shown  in  Figure  3  as  well.  The  important  trend  made  visible 
is  this  investigation  as  shown  in  Figure  3,  is  that  k^^  (and  hence 
decreases  with  increasing  nondimensionalized  coupling  L2 , 
which  perhaps  may  be  a  more  effective  way  to  actually  measure  and 
compare  elastic  couplings,  and  ^2S'  example,  the  results 

from  MIT^'^  which  were  computed  for  some  representative  configura¬ 
tions,  in  dimensionalized  form  seem  to  represent  systems  with  a 
fairly  significant  variation  in  coupling  D2g  or  (depending  on 

the  coordinate  system) .  However  when  nondimensionalized,  the  re¬ 
sults  as  shown  in  Figure  3  seem  to  show  little  variation  in  coup¬ 
ling.  In  fact  they  appear  to  be  so  close  to  the  zero  coupling 
case  (X  /A  =  1)  or  isotropic  (or  metal)  case,  that  k^^  for  a  metal 
or  isotropic  wing  should  be  a  good  approximation  (if  it  was  nec¬ 
essary  to  make  an  approximation) .  The  low  value  for  the  effective 
coupling  was  also  evident  in  the  nondimensionalized  results  from 
jip3 , 4  the  bending  freguency  hardly  _yaried  with  material 

changes.  The  question  that  could  be  asked  is  therefore  "Do  all 
possible  composite  wing  configurations  result  in  very  low  effective 
nondimensionalized  coupling,  (A^/A^=>1)7"  If  the  answer  is  yes  , 
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then  it  may  be  proposed  that  for  A  >.5,  k,  (A  ,A  /A  )  may  be  ap- 
proximated  as  k^(A^,l),  which  is  also  the  isotropic  or  metal 
value.  For  this  case  it  is  seen  that  the  computation  of  the  na¬ 
tural  frequencies  of  composite  aircraft  wing  having  aeroelastic 
oscillations,  merely  requires  the  computation  of  (D22/P)  * 

given  wing  half  span  I  ,  since  k  (A  ,A  /A  )  is  approximately 
equal  to  k  (A  ,1)  which  is  approximately  equal  to  a  constant 
(3.5)  for  n  =  1.  This  result  should  make  frequency  computations 
for  the  bending  mode  significantly  easier. 

The  problem  with  an  affirmative  answer  which  may  likely  be 

a  "practical"  answer,  to  the  question  posed  above,  is  that  there 

doesn't  seem  to  be  a  theoretical  or  rigorous  analytical  reason 

(to  the  best  of  the  author^'  knowledge)  why  X  /A  must  always  be 

c  c 

approximately  1.  Therefore,  if  on  the  other  hand,  the  answer  to 

our  question  is  negative,  then  the  following  observations  may  be 

mades  (a)  significant  variation  in  k^  is  possible  with  variations 

in  effective  nondimensionalized  coupling  (A  /A  ) .  In  fact  it  can 

c  c 

be  seen  from  Figure  3  that  if  X^/A^  approaches  zero,  k^^  (and  hence 

tOj^)  approaches  zero.  (b)  The  values  of  kj^  vary  significantly  with 

A  for  low  A  but  approach  asymptotic  values  for  large  A  . 
c  C 

(c)  The  highest  values  of  k^^  is  for  isotropic  (metals)  or  quasi¬ 
isotropic  configurations.  (d)  For  large  values  of  A^,  there  ap¬ 
pears  to  be  a  simple  approximate  (hopefully  linear  relationship) 
between  k^^  and  X^/A^  (or  a  measure  of  coupling).,  (e)  For  very 
large  coupling  (X^/A^->0)  ,  kj^  approaches  zero,  which  may  provide 
the  ingredient  necessary  for  coupling  between  the  elastic  and  ri¬ 
gid  motions. 
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Perhaps  a  number  of  implications  of  some  of  these  observations 
should  be  examined:  Observation  (c)  seems  to  imply  that  the  high¬ 
est  first  frequency  would  correspond  to  isotropic  or  quasi-isotrop¬ 
ic  configurations  if  (D22/P)  and  are  the  same.  It  is  known  how¬ 
ever  the  metals  have  lower  values  of  (D22/P)  than  composites.  It 
therefore  means  that  quasi-isotropic  or  orthotropic  configurations 
are  desirable  for  such  a  design  goal.  Observation  (e)  would  seem 
to  imply  that  if  a  designer,  interested  in  tailoring  the  wing  fre¬ 
quencies,  arbitrarily  introduces  large  effective  nondimensionalized 
coupling  »  then  (and  hence  would  approach  zero. 

This  may  result  in  coupling  between  flexible  and  rigid  motions 
which  may  or  may  not  lead  to  instabilities.  Could  this  have  hap¬ 
pened  in  the  case  of  the  X-29  Forward  Swept  Composite  Wing  Aircraft 
for  which  one  of  the  primary  modes  of  instability  results  from  the 
coupling  between  flexible  and  rigid  body  motions?  In  other  words 
was  "too  much"  coupling  (effective)  inadvertently  built  into  the 
wing  during  the  design  process?  If  that  is  the  case,  is  there  an 
alternative,  equivalent  design  without  any  penalties  (weight  or 
otherwise)  that  could  have  been  explored?  Although  the  answers  to 
these  questions  can,  strictly  speaking,  only  be  possible  after  car¬ 
rying  the  necessary  aeroelastic  analysis  in  which  unsteady  aerody¬ 
namic  forces  are  considered,  it  appears  from  Figure  3  that  a  rough 
idea  of  the  final  picture  may  be  obtained  from  the  natural  fre¬ 
quency  analysis.  Afterall,  it  is  a  common  belief  that  the  phenom¬ 
ena  that  actually  lead  to  aeroelastic  instabilities  are  linked  to 
damping  and  coupling. 
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Earlier  in  this  paper  it  was  mentioned  that  previous  studies 
by  other  investigators  have  found  what  appeared  to  be  some  kind 
of  modal  transformations  as  ply  orientation  was  changed  in  a  de¬ 
sign  process  for  a  composite  wing.  While  variation  in  ply  orien¬ 
tation  may  change  several  directional  stiffness  parameters  for 
the  wing,  the  coupling  stiffness  parameter,  L2 ,  may  be  singled 
out  as  a  significant  design  parameter,  because  it  may  vary  con¬ 
siderably  (orthotropic  configurations  have  zero  values  while  it 
may  be  fairly  significant  in  other  configurations) .  Furthermore 
it  should  be  remembered  that  the  main  reason  for  ply  orientation 
variation  is  for  'tailoring',  which  is  believed  to  be  primarily 
tied  to  couplings  and  ^2^)  •  The  absence  (or  the  presence) 

of  these  couplings  is  basically  what  differentiates  orthotropic 
configurations  from  anisotropic  configurations.  From  these  ob¬ 
servations,  and  the  fact  that  the  entity  that  ties  the  bending 
and  torsional  equations  is  the  coupling,  it  became  clear  that 
the  role  of  coupling  in  modal  transformations  should  be  signifi¬ 
cant. 

In  order  to  see  the  role  of  coupling  therefore  in  this  stu¬ 
dy,  the  modal  assumptions  for  the  coupled  problem  were  made  sim¬ 
ilar  to  those  normally  made  for  the  uncoupled  problem  (e.g.,  the 
frequency  was  assumed  to  be  real)  so  as  to  provide  an  opportun¬ 
ity  to  compare,  contrast  and  discern  the  final  results  easily. 
When  this  was  done  and  the  eigen-problem  was  formulated  result¬ 
ing  in  a  complex  transcendental  expression  from  which  the  eigen¬ 
value  are  to  be  extracted,  a  careful  examination  began. 

A  significant  difference  between  the  coupled  and  uncoupled 


problems  was  that  (as  shown  in  Figure  2)  the  transcendental  ex¬ 
pression  from  which  the  eigenvalues  are  extracted  is  complex  for 
the  coupled  problem  while  it  is  real  for  uncoupled  problems. 

The  complex  nature  of  this  transcendental  expression  basically 
reflects  the  fact  that  the  bending  and  twisting  oscillations  are 
generally  out  of  phase.  Therefore  the  resultant  coupled  frequen¬ 
cy  that  represents  both  the  bending  and  twisting  oscillation  may 
be  viewed  as  some  kind  of  vector  representation  of  the  individual 
contributions.  In  order  to  formulate  some  explanations  for  the 
phenomenon  of  modal  transformation  in  coupled  (conservative)  sys¬ 
tem,  it  may  be  necessary  to  compare  and  contrast  coupled  systems 
with  damped  systems.  Damped  systems,  by  definition  are  noncon¬ 
servative,  i.e.,  the  system  experiences  a  net  loss  or  gain  in  en¬ 
ergy.  It  is  well  known  that  in  a  damped  system  the  transcenden¬ 
tal  expression  for  extracting  the  eigenvalues  is  complex,  again, 
due  to  the  phase  angle  that  exists  between  the  damping  force  and 
the  conservative  forces  in  the  system.  It  is  also  known  that 
some  damping  (desirable  types)  would  tend  to  reduce  the  oscilla¬ 
tion  of  the  system  (the  non-desirable  type  tend  to  make  the  oscil¬ 
lations  diverge) .  Therefore  since  damping  is  linked  to  some  ener¬ 
gy  transfer  which  in  turn  tend  to  lead  to  a  change  in  the  oscilla¬ 
tion  frequencies,  it  was  thought  that  the  complex  nature  which  is 
common  to  the  coupled  and  d2unped  system  determinants  (transcenden¬ 
tal  functions)  from  which  the  eigenvalues  are  extracted  may  be  a 
similarity  that  may  provide  some  explanations  to  the  modal  trans¬ 
formations  in  coupled  systems.  Using  the  similarity  argument,  the 
coupled  system  which  for  the  present  problem,  is  conservative,  may 


71 


be  viewed  as  having  some  conservative  inter-modal  energy  trans¬ 
fer  within  the  system  when  the  coupling  is  changed,  resulting 
in  steady  state  changes  or  transformations  of  the  modal  energy 
content  of  a  coupled  mode  compared  to  the  uncoupled  case.  It 
may  be  worthwhile  to  point  out  that  some  results  recently  ob~ 

Q 

tamed  at  Purdue  University  and  communicated  to  the  authors 
seem  to  strongly  support  this  hypothesis. 


5.3.10 

Empirical  Relations 

A  careful  study  of  Figure  3  has  led  the  authors  to  propose 
the  following  closed  form  asymptotic  relationship  that  may  be 
useful  for  some  preliminary  design  consideration: 

ki  =  3.5  ^ 

Equation  28  was  derived  from  Figure  4.  Equation  28  as  well  as 

equation  27  show  that  the  first  coupled  frequency  decreases  with 

increasing  coupling,  a  trend  that  seems  to  be  supported  by  new 

8  3 

results  from  Purdue  University  and  the  data  from  MIT.  In  ref¬ 
erence  3  for  example,  the  first  nondiroensionalized  frequency  com¬ 
puted  by  Raleigh-Ritz  (in  which  coupling  is  zero)  had  a  value  of 
3.52,  which  is  consistently  higher  than  those  computed  by  finite 
element  method  in  which  coupling  is  finite  (not  equal  to  zero) . 
Equations  27  and  29  which  are  closed  form  (generally  rare  for  an- 
istropic  systems)  should  be  easy  to  use. 
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Before  this  discussion  is  concluded,  it  is  probably  necess¬ 
ary  to  explain  why  only  the  results  of  the  first  mode  are  shown 
in  this  paper.  First  of  all  it  should  be  pointed  out  that  some 
second  mode  data  have  been  generated  but  are  still  being  studied 
very  critically  to  understand  the  general  trends.  It  may  also 
be  pointed  out  that  the  extraction  of  the  eigenvalues  is  a  little 
challenging  since  some  care  is  needed  in  ordering  the  roots  of 
the  operator  equations. 

5.3.11 

Concluding  Remarks 

This  paper  has  attempted  to  present  exact  closed  form  solu¬ 
tions  to  the  coupled  bending-torison  vibration  problem  for  a 
simplified  model  of  composite  aircraft  wings  with  war  ping  con¬ 
straint.  Increasing  the  coupling  was  found  to  decrease  the  first 
coupled  frequency.  A  comparison  between  the  coupled  problem  and 
a  sample  damped  problem  .led  the  authors  to  propose  some  explana¬ 
tion  to  the  "modal  transformation"  phenomenon  found  by  earlier  in¬ 
vestigators.  Some  simplified  closed  form  expressions  are  provi¬ 
ded  for  the  first  coupled  frequencies  which  may  be  useful  for  fast 
applications . 
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5.4  INVESTIGATION  OF  THE  BASIC  MECHANISM  OF  AEROELASTIC  INSTABILITY 


FOR  COMPOSITE  AIRCRAFT  WINGS 

It  is  well  known  in  aeroelastici ty  that  the  mechanism  of  instabil¬ 
ity  is  governed  by  the  interaction  of  aerodynamic  forcing  functions 
and  the  vibrational  modes  of  the  structure.  The  aerodynamic  forces 
generally  provides  ingredients  such  as  damping  and  coupling  to  such 
an  interaction.  As  a  result,  the  structural  modes  may  loose  damp¬ 
ing  and/or  stiffness  which  may  result  in  some  form  of  coalescence 
of  two  or  more  modes  leading  eventually  to  instability  which  may  or 
may  not  be  catastrophic.  Physically  therefore  it  is  sensible  to  say 
that  some  of  the  main  ingredients  of  aeroelastic  instability  are 
damping  and  coupling. 

Crisp^  has  shown  that  the  necessary  condition  for  instability  may  be 
studied  by  investigating  the  matrix  form  of  the  aeroelastic  equa¬ 
tions  of  motion.  Thus 

[A] {q^}  +  [B]  {q. }  +  [C]  {q^}  =  0  (1) 

where  q^  are  the  mode  shapes  and  A,B,C,  are  the  mass,  damping  and 
stiffness  matrices,  respectively.  He  considered 

T  =  1/2  [q^]  [A]  {q^} 

F  =  1/2  [q.]  [B]  {q.}  (2) 

U  =  1/2  [q^]  IC]  {q^} 

where  i  and  U  are  some  generalized  kinetic  and  potential  energies, 

2 

respectively.  F  is  basically  the  Rayleigh's  dissipation  function  . 
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(A]  ,  [B]  and  IC]'  are  symmetrical  for  elementary  metal  mechanical 

systems  but  the  addition  of  the  aerodynamic  and  dampings  spoils 
this  symmetry  for  metal  or  isotropic  systems  for  the  last  two.  In 
the  case  of  composite  anisotropic  system,  aerodynamics  is  not  even 
needed  for  asymmetric  [C] ,  since  such  asymmetry  can  be  provided  by 
elastic  couplings.  Crisp  has  further  shown  in  the  case  with  non- 
synunetric  [B]  and  [C]  they  may  be  decomposed  into  symmetric  and 


skew-  symmetric 

portions 

or 

[B]  = 

[B^l  + 

IB2] 

(3) 

[C]  = 

[C^]  + 

IC2J 

where  IB^]  ,  (C^l  are  the  symmetric  portions  and  tB2]  ,  IC2I  are  the 
skew- symmetric  portions.  By  considering  the  total  energy  for  the 
system  t  +  U  =  E,  Crisp  obtained  the  following  expression  for  work 
done  on  this  system, 

~  =-  Iq^l  IB^]  {c^.}  -  Iq.)  IC2]  {q^}  (4) 

It  is  then  seen  that  without  damping  (  [Bj  =  0)  ,  a  composite  wing 
can  exhibit  a  "damping- like"  behavior  as  was  shewn  in  section  5.4 
probably  for  the  first  time  in  aeroelasticity .  This  investigation 
is  being  continued  in  order  to  search  for  the  explanation  for  aero- 
elastic  tailoring.  A  preliminary  plot  of  the  second  mode  is  pre¬ 
sented  in  Figure  1.  Figure  2  also  shows  that  modal  coalescence 
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is  also  a  reality.  Further  studies  have  also  shown  that  such  a 
coalescence  could  lead  to  instability.  We  are  currently  exam¬ 
ining  the  necessary  physical  ingredients  for  such  an  instabil¬ 
ity  and  its  implications.  Figure  3  shows  a  typical  flutter  in¬ 
stability  evolution  process.  Figures  2  and  3  therefore  display 
some  interesting  resemblance.  We  are  in  the  process  of  discern¬ 
ing  the  striking  resemblance. 
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Nomenclature 

^  piiysical  and  affine  space 
coordinates,  respectively 

>  chordwise  Integrals 

■  differential  aerodynamic  pressure  dis¬ 
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=  displacement 
=  wing  box  depth 

=  affine  space  bending  and  torsional  dis¬ 
placement,  respectively 
=  affine  space  half-chord  and  chord,  res¬ 
pectively 

=  affine  space  half-span  for  the  wing 
=  parameter  that  measures  the  location  of 
the  elastic  axis  relative  to  mid-chord 
=  affine  space  mass  per  unit  span 
=  vibration  frequency 


Introduction 

The  performance  of  modern  supermaneuverable  aircraft  can  be  made 
to  benefit  a  great  deal  from,  significant  advances  in  materials 
technology  and  the  availability  of  more  accurate  aerodynamic  pre¬ 
diction  capabilities.  Supermaneuverability  as  a  design  goal  in¬ 
variably  calls  for  an  optimization  of  the  design  parameters.  Op¬ 
timization  may  be  partially  accomplished  for  example,  by  using 
composite  materials  to  minimize  weight.  Indeed,  it  has  been  known 
that  these  composite  materials  can  be  tailored  to  resolve  the  dy¬ 
namic  or  static  instability  problems  of  these  types  of  aircraft. 
The  concept  is  referred  to  as  aeroelastic  tailoring. 


While  aeroelastic  tailoring  has  tremendous  advantages  in  the  design 
of  an  aircraft,  the  analysis  which  provides  the  basis  for  the  aero¬ 
elastic  tailoring  itself  is  generally  very  involved.  This  is  ra¬ 
ther  unfortunate  since  a  good  fundamental  physical  insight  of  the 
tailoring  mechanism  is  required  for  accurate  and  reliable  results. 
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In  this  paper  an  attempt  is  made  to  look  at  some  dynamic  theories 
that  can  be  used  to  understand  the  aeroelastic  tailoring  mechan¬ 
ism.  Specifically,  the  accuracy  of  the  St.  Venant  torsion  theory 
which  is  relatively  simple  and  frequently  used  in  aeroelastic 
analysis  is  examined  with  particular  reference  to  the  effects  of 
the  wing's  aspect  ratio  as  well  as  other  design  parameters. 

12  3 

Although  earlier  studies  '  '  have  indicated  that  the  St.  Venant's 
torsion  theory  is  reasonably  accurate  except  for  aircraft  wings 
with  fairly  low  aspect  ratios,  the  theory  supporting  that  conclu¬ 
sion  was  based  on  the  assumption  that  the  wing  is  constructed  of 
isotropic  materials.  Basically,  the  St.  Venant's  torsion  theory 
assumes  that  the  rate  of  change  of  the  wing's  twist  angle  with 
respect  to  the  spanwise  axis  is  constant  { for 'constant  stiffness 
and  torque) .  This  assumption  is  hardly  accurate  particularly  for 
modern  aircraft  construction  in  which  different  construction  mate¬ 
rials  are  employed  and  the  aerodynamic  loads  vary  significantly 
along  the  wing's  span^  However,  References  1-3  have  shown  that 
(in  spite  of  such  an  inaccurate  assumption)  the  main  parameter 
that  determines  the  accuracy  of  the  St,  Venant's  theory  is  the 
wing's  aspect  ratio.  Thus,  it  was  determined  that  the  theory  is 
fairly  accurate  for  moderate  to  high  aspect  ratio  wings  construc¬ 
ted  of  isotropic  materials.  In  recent  studies^ however,  it  has 
been  shown  that  for  wings  constructed  of  orthotropic  composite 
materials,  the  conclusions  of  References  1-3  need  to  be  modified. 
Rather  than  using  the  geometric  aspect  ratio  of  the  wing  to  deter¬ 
mine  the  accuracy  of  St.  Venant's  twist  theory,  it  was  suggested 
that  a  generic  stiffness  ratio,  as  well  as  an  effective  aspect 
ratio  which  considers  the  wing's  geometry  and  the  ratio  of  the 


principal  directional  stiffness,  should  be  considered  in  estab¬ 
lishing  the  accuracy  of  St.  Venant's  theory. 


The  present  paper  is  basically  an  extension  of  the  studies  that 
were  begun  in  References  4  and  5.  In  thi-  study  the  task  was  to 
examine  the  role  of  coupling  (elastic  coupling)  on  the  accuracy 
of  St.  Venant's  theory  applied  to  static  problems.  It  was  dis¬ 
covered  that  coupling  plays  a  very  significant  role  on  the  ac¬ 
curacy  of  St.  Venant's  twist  theory. 


5.5.3 


Formulation 


Consider  an  aircraft  wing  fabricated  of  composite  materials  and 
mathematically  idealized  as  a  canti-levered  plate  subjected  to 
forces  and  moments.  It  can  be  shown  that  the  equations  of  motion 
for  such  a  model  can  be  described  as  follows. 
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D. ,  are  the  elastic  constants,  p,  is  the  material  density,  Ap  is 
1  j 

the  differential  pressure  distribution,  vr  is  the  displacement,  t 
is  the  time,  and  h  is  the  wing  box  depth. 

5.5.4 

Evolution  of  Warping  Parameters 

The  evolution  of  the  warping  parameter  with  which  to  study  the 
aeroelastic  warping  constraint  phenomenon  for  wings  fabricated 
of  composite  materials  is  a  process  that  depends  on  the  sophisti¬ 
cation  of  the  wing's  mathematical  model;  whether  coupling  ef¬ 
fects  are  included,  whether  the  wing's  chordwise  curvatures  are 
included  and  so  on.  Therefore,  any  warping  parameter  is  as  good 
as  the  corresponding  wing's  displacements  assumptions.  However, 
the  virtual  work  equation  makes  it  possible  for  the  analyst  to 
determine  its  effective  independent  variables  even  before  the  dis- 
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placement  assumptions  are  made.  By  non-dimensionalizing  the 
spanwise  space  variable  in  Equation  (1) ,  depending  on  whether  one 
is  interested  in  the  static,  dynamic,  coupled  or  uncoupled  dis¬ 
placements,  one  of  the  following  warping  parameters  may  be  useful. 


where 

D*  =  D*{l-e)  (7) 

(£.^/c^)  is  defined  as  the  wing's  effective  aspect  ratio  and 
and  L  are  the  generalized  stiffness  and  coupling  ratios  respective¬ 
ly  (defined  in  earlier  work  such  as  References  5  and  6) , 

Equations  (5)  and  (6)  represent  the  appropriate  warping  parameter 
for  dynamic  deformation,  static  displacement  with  elastic  cross¬ 
coupling. 

It  was  discovered  in  this  study  that  evolving  the  warping  param¬ 
eter  in  a  manner  shown  in  Equation  (1)  thru  (3) ,  should  enable  one 
to  effectively  investigate  the  effects  of  warping  on  the  composite 
wing's  dynamics  (or  the  accuracy  of  St.  Venant's  theory.)  From  the 

ft 

Icunination  theory  for  composites  it  is  known  that  while  and 
(i  /c  )  are  always  positive,  L  can  be  positive  or  negative. 

However,  from  Equations  (1)  and  (3) ,  it  is  clear  that  whether  a 
composite  wing  has  positive  or  negative  coupling,  the  warping  ef¬ 
fect  (in  terms  of  X^)  is  unchanged. 
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5.5.5 


Computations 

By  using  the  evolved  warping  parameters  defined  in  Equations  7  and 
8  and  appropriate  boundary  conditions,  the  boundary  value  problems 
associated  with  Equation  (4)  are  solved  in  a  closed-form  manner  to 
determine  the  wing's  static  twist. 

The  wing  loading  conditions  considered  in  this  analysis  are  as 
follows:  (a)  steady  state  distributed  twist  loads  and  (b).  steady 

state  concentrated  twist  oads. 

5.5.6 

a.  Steady  distributed  twist  loads 

For  a  wing  with  a  constant  uniformly  distributed  spanwise 
twisting  moment,  f^  resulting  from  a  steady  state  coupled 
bending-torsion  displacements,  the  exact  closed  form  solutions 
for  the  mode  shape  satisfying  the  appropriate  boundary  con¬ 
ditions  is  given  by 
,2 
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Equation  8  is  therefore  a  closed  form  coupled  twist  distribution 
for  a  composite  wing  with  the  warping  effects  accounted  for. 
where 

fo  - 

When  equation  8  is  evaluated  at  the  wing  tip  and  compared  to  an 
equivalent  expression  predicted  by  St.  Venant's  theory,  the  fol¬ 
lowing  expression  is  obtained. 
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where  (1)  is  the  wing  tip  twist  given. by  equation  15  while 
(1)  is  the  wing  tip  twist  given  by  the  St.  Venant  torsion 

theory.  A  plot  of  equation  10  is  shown  in  Fiqure  1. 

5.5.7 

b.  Steady  state  concentrated  tip  twist  loads 

If  the  wing  is  under  the  influence  of  a  concentrated 
twisting  moment,  at  the  tip  as  a  result  of  a  steady 
state  coupled  bending  torsion  displacement,  the  exact 
closed  form  twist  distribution  that  satisfies  these 
equations  of  motion  and  their  associated  boundary 


conditions  is  given  by 


.  o  o  — 

o  c  L 


(coshA^y^  -1) 


sinh4A^y  tanh4A 
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j 

When  the  twist  distribution  given  by  equation  11  is  evaluated  at 
the  wing  tip  and  compared  to  its  counterpart  predicted  by  the 
St.  Venant 's  torsion  theory  the  following  expression  is  obtained, 


“o^^^  St.V 


tanh4 A 


=  1- 


It  should  be  noted  that  the  ratio  given  by  equation  12  was  plot¬ 
ted  for  the  real  values  of  A^  in  references  5  and  6  and  was  shown 
to  represent  conditions  where  any  errors  resulting  from  using  St. 
Venant' s  torsion  theory  are  conservative  (over-design  rather  than 
under-design) .  In  this  analysis  equation  12  is  examined  when 
is  imaginary,  which  is  possible  if  is  very  large.  Under 
such  circumstances,  equation  12  becomes 
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Figure  2  depicts  the  conditions  given  by  equation  13.  It  is 

therefore  seen  from  the  figure  that  there  are  certain  ranges  of 

for  which  nonconservative  errors  are  possible  by  using  the 

St.  Venant's  twist  theory. 

5.5.8 

Results  and  Conclusions 

The  results  are  shown  in  Figures  1  and  2.  Figure  1  shows  a  com¬ 
parison  of  the  static  wing  tip  twist  obtained  in  the  present 
study  and  that  obtained  via  St.  Venant's  twist  theory  in  the  pres¬ 
ence  of  statically  distributed  forces  and  low  to  moderate 
coupling.  Figures  2  shows  the  trend  for  concentrated  forces  and 
substantial  coupling.  In  Figure  1  it  is  seen 

that  the  presence  of  coupling  makes  the  errors  of  St.  Venant's 
theory  worse.  This  seems  to  suggest  that  the  more  sophisticated 
theory  is  more  important  for  wings  with  coupling  (e.g.,  wings 
aeroelastically-tailored  using  elastic  cross-coupling) . 

Figure  2  also  shows  that  nonconservative  errors  (] — pr-r-  j>l)  are 

®o '  ' 

are  possible.  St.V 

Using  Figures  1  and  2,  the  following  conclusions  can  be  summar¬ 
ized:  (i)  ignoring  warping  arbitrarily  using  St.  Venant's  theory 

could  result  in  very  significant  errors  (as  high  as  over  80%  er¬ 
rors)  in  analytical  results  for  composite  aircraft  wings,  (ii) 
warping  is  more  important  (St.  Venant’s  theory  is  less  accurate) 
for  wings  with  coupling,  (iii)  St.  Venant's  theory  (which  has  al- 
ways  been  shown  to  be  conservative'  '  '),  can  be  non-conservative 
or  St.  Venant's  approximation  can  lead  to  an  unsafe  design  error 
(under  design  rather  than  over  design  from  a  stability  point  of 


view)  . 
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Figures 


Figure  1 


Figure  2 


Wing  Tip  Twist  Ratios  for  Simple  and  More  Involved 
Theory  (distributed  load  and  low  to  moderate  coupling) 


Wing  Tip  Twist  Ratios  Comparing  Simple  and  More  In¬ 
volved  Theory  (with  substantial  coupling) 
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Closed-Form  Solutions  for  Nonlinear  Quasi-Unsteady 
Transonic  Aerodynamics 
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The  existtiice  of  exact  dosed-fonn  solotiom  Tor  nonlinear  unsteady  aerodynamics  is  estabUshed.  The  full 
nonlinear  unsteady  velocity  potential  equations  for  an  airfoil  are  considered.  Evidence  Indicating  why  the  traditional 
hodograpb  approach  is  ineffective  for  solving  these  equations  is  provided.  Therefore,  a  suitable  mapping  scheme  is 
employed  in  transforming  these  full  nonlinear  equations  into  the  hodograph  plane.  A  close  examinatioo  of  the 
resulting  hodograph  equations  reveals  closed-form  solntioos  can  be  obtained  for  the  nonlinear  nnsteady  aerodynamic 
characteristics  of  an  airfoil  in  potential  flow.  The  sbockless  transonic  results  presented  in  this  hvisdd  analysis  show 
trends  that  are  in  agreement  with  the  results  of  previous  investigators  and  available  expeiimental  data.  “Dips”  were 
observed  In  the  pressure  distrilmtloos  as  the  freestieam  Mach  number  is  varied.  It  appears  that  there  are  Hnite 
optimum  redneed  frequencies  for  the  pressure  distributions.  This  result  might  suggest  a  solntioo  to  the  “transonic 
dip”  problem.  Perhaps  an  important  practical  consequence  of  this  study  b  the  possibility  of  employing  thb 
approach  to  solve  an  inverse  problem  of  designing  an  airfoil  section  with  given  or  desired  aerodynamic 
characterfadcs.  Desirable  candidates  for  such  a  design  procedure  wonld  include  superciitical  oscOlatiiig  shock-free 
or  “transonic  dipless"  airfoil  sections.  Such  airfoils,  therefore,  could  be  designed  to  meet  both  the  performance  and 
stability  criteria  simnltaneously. 


Nomenclature 


b,c 


(x,?):  (u.tv) 
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q.P 
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Q^Q 
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p,m 


f.y 
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=  airfoil  chord  length  and  velocity  of  sound, 
respectively 

=  Cartesian  and  hodograph  coordinates, 
respectively 

=  pressure  and  lift  coefRcients,  respectively 
=  reduced  frequency  and  Mach  number, 
respectively 

=  resultant  flow  velocity  and  pressure, 
respectively 

=  Jacobian  of  the  hodograph  transformation 
=  hypergeometric  function  and  arbitrary 
constant,  respectively 
=  velocity  potential  and  stream-function 
quantity,  respectively 
=  integer  and  square  root  of  minus  one, 
respectively 

=  air  density  and  oscillation  frequency, 
respectively 

=  time  and  ratio  of  gas  specific  heats, 
respectively 

=  velocity  potential  and  vector  differential 
operator,  respectively 

=  transformed  velocity  potential  and  stream 
function,  respectively 
=  affine  space  or  nondimensionalized 
quantities 

=  quantities  at  infinity 
=  angle  inclined  by  velocity  vector  and 
positive  X  axis 

=  dimensionless  velocity  variable 
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Introduction 

The  hodograph  transformation  has  been  established  as  a 
very  vital  tool  with  which  to  analyze  two-dimensional 
nonlinear  steady  transonic  flow  problems  comprehensively. 
The  significance  of  the  hodograph  representation  results  from 
the  fact  that  whereas  the  linearized  equations  in  the  physical 
plane  fail  to  explain  certain  observed  transonic  flow  phenom¬ 
ena.  the  nonlinear  equations,  which  more  accurately  dwrite 
the  flow,  do  not  seem  to  have  simple  closed-form  solutions  in 
the  physical  plane.  With  the  help  of  suitable  mapping  func¬ 
tions,  these  nonlinear  equations  could  be  transformed  into  the 
hodograph  plane,'  '*  where  they  become  linear,  resulting  in  a 
possibility  of  obtaining  closed-form  solutions  for  them.  Al¬ 
though  the  hodograph  approach  could  introduce  difficulties 
associated  with  involved  boundary  conditions  resulting  from 
some  practical  problems,  many  investigators  believe  that  its 
advantages,  particularly  in  terms  of  physical  insights,  far 
outweigh  its  disadvantages. 

Outstanding  contributions  from  various  investigators  have 
been  responsible  for  the  development  of  the  hodograph 
method  for  solving  nonlinear  potential  flow  problems.  The 
early  work  of  Molenbrock'  and  Chaplygin*  was  complemented 
by  the  later  efforts  of  investigators  like  LighthilP  and  Guder- 
leg  *  Consequently,  researchers  like  Nieuwland,’  Bauer, 
Garabcdian.  and  Korn,*  and  Boersteel*  were  able  to  establish 
the  hodograph  approach  as  an  effective  design  tool  for 
efficient  airfoils  like  the  supercritical  shockless  sections.  The 
basic  idea  here  is  to  suppress  the  boundary-layer  separation 
by  "pushing”  the  shock  waves  on  the  wing  toward  its  trailing 
edge  and  eventually  diluting  (or  weakening)  them  as  much  as 
possible.  Sobieezky’  ’  and  his  collaborators  have  also  pre¬ 
sented  interesting  results  more  recently. 

The  hodograph  transformation  has  been  known  and  used 
for  over  a  century.  Curiously,  however,  evidence  from  a 
literature  search  seems  to  indicate  that  its  use  has  largely  been 
restricted  to  the  analysis  of  the  nonlinear,  steady,  two-dimen¬ 
sional  flow  problems  like  steady  transonic  flow. 

The  reasons  for  this  restriction  or  why  there  has  not  been 
an  extension  of  this  approach  to  unsteady  two-dimensional  or 
three-dimensional  flow  problems  also  appear  to  be  absent  in 
the  literature  as  well.  The  implication,  therefore,  seems  to  be 
that  such  an  approach  can  only  be  used  to  solve  steady-state 
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airfoil  problems  like  performance  and  not  dynamic  instability 
problems  (which  require  the  unsteady  flow  solutions). 

In  transonic  aeroelasticity,  phenomena  like  the  "transonic 
dip"  require  nonlinear  unsteady  aerodynamics  for  their 
understanding.  The  need  to  operate  modem  aircraft  in  the 
transonic  region,  among  other  things,  has  been  responsible  for 
the  recent  tremendous  interest  in  transonic  flow  problems. 
Aeroelasticians  and  computational  fluid  dynamicists  have 
been  investigating  these  problems  with  a  good  degree  of 
success;  Refs.  15-24  are  just  a  few  examples.  As  a  result, 
numerous  computer  codes  are  now  available  for  computing 
transonic  aerodynamics,  many  of  them  using  various  approx¬ 
imations  of  the  full  nonlinear  potential  equations.  In  spite  of 
this  progress,  however,  recent  publications  and  presentations, 
like  Ref  24,  seem  to  indicate  that  a  great  deal  still  remains  to 
be  done  to  thoroughly  understand  aeroelastic  phenomena  like 
the  transonic  dip.  The  author,  inspired  by  a  recent  experi¬ 
ence,”  shares  the  view  held  by  some  investigators  that  a  more 
fundamental  approach  could  provide  some  of  the  necessary 
physical  insights  for  improving  our  understanding  of  these 
phenomena.  Such  insights  may  not  be  readily  extractable 
from  an  analysis  in  which  big  computer  codes  are  used. 

This  paper  presents  the  results  of  a  preliminary  investiga¬ 
tion  of  the  unsteady  transonic  aerodynamics  in  which  the  full 
nonlinear  unsteady  two-dimensional  velocity  potential  flow 
equations  are  employed.  Evidence  indicating  why  the  tradi¬ 
tional  hodograph  approach  is  not  effective  for  solving  these 
equations  is  presented.  With  the  help  of  a  certain  mapping 
scheme,  the  nonlinear  equations  are  transformed  into  the 
hodograph  plane.  A  close  examination  of  the  transformed 
equations  reveals  that  if  7,  the  Jacobian  of  the  transformation, 
is  prescribed,  ab  initio,  the  hodograph  velocity  potential 
satisfies  a  linear  second-order  partial  differential  equation. 
The  significance  of  this  result  includes  the  fact  that  the 
following  results  appear  to  be  possible  for  the  first  time: 

1)  Exact  closed-form  solutions  for  the  nonlinear  unsteady 
velocity  potential  can  be  obtained. 

2)  The  solution  of  an  inverse  problem  of  designing  an 
airfoil  section  with  given  or  desired  unsteady  aerodynamic 
characteristics  can  be  attempted. 

These  results,  therefore,  imply  that  it  is  possible  to  design  an 
airfoil  that  can  meet  both  the  performance  and  dynamic 
stability  criteria  simultaneously. 

The  results  shown  in  this  paper  are  obtained  by  piecing  the 
fundamental  solutions  in  a  manner  similar  to  Nieuwland's 
approach.’  From  these  results  it  appears  that  there  are  "dips" 
in  the  pressure  distributions  as  the  freestream  Mach  number  is 
varied  in  the  transonic  region,  a  phenomenon  that  also  has 
been  established  (Ref  24,  for  example)  for  the  aeroelastic 
stability  characteristics.  It  also  is  observed  that  there  appear 
to  be  finite  optimum  reduced  frequencies  for  the  pressure 
distributions.  This,  therefore,  appears  to  agree  with  the  trends 
established  by  Marble”  for  a  quasi-one-dimensional  flow. 

Equations  of  Motion 

By  considering  the  concept  of  "control  volume"  for  a  fluid 
flow,  it  can  be  proved  from  the  first  principles  that  the 
physical  principle  of  continuity  of  mass  demands  that  the 
following  equation  be  satisfied  for  any  arbitrary  volume: 

(^+?  v)p-(-p(V  »)=0  (1) 

where  f  is  the  fluid-flow  velocity  vector,  p  is  the  fluid  density, 
and  V  is  a  differential  vector  operator. 

Using  a  similar  approach,  the  physical  principle  of  the 
conservation  of  momentum  can  be  represented  by  the  follow¬ 
ing  equation: 

+(?  V)j|f)v ?)  = -VF  (2) 


By  using  Leibnitz's  rule  and  Kelvin’s  theorem  for  an  irrota- 
tional  flow  as  well  as  the  assumption  that  a  velocity  potential 
<f)  exists,  such  that 


9  =  V0  (3) 

in  light  of  Eq.  (1),  Eq.  (2)  can  be  integrated  to  obtain  the 
following: 

#  V0  f"  dP 

-57+ - z - +  — =-:r-  (4) 

2  p  2 

where  ( ),,  is  a  reference  point  quantity  (e.g.,  infinity).  Equa¬ 
tion  (4)  is  basically  the  l^rnoulli- Kelvin  equation. 

By  considering  isentropic  flow  with  the  following  relation¬ 
ship  between  the  pressure  P  and  density 


pp  ’  =  const 
and  defining  the  quantity  c  given  by 


(5) 


(6) 


called  the  speed  of  sound,  the  Leibnitz  rule  can  be  employed 
to  derive  the  following  velocity  potential  equation  from  Eqs. 
(I)  and  (4): 


- 1^^  (V4>  +  ^  +  Vm  ■  V^)  J  =  0  (7) 

When  Eq.  (7)  is  expanded  and  the  vector  calculus  is  carried 
out,  the  following  nonlinear  velocity  potential  equation  in  the 
physical  Cartesian  coordinates  is  obtained: 


(c*-u*)  .^-2uh’ 


SxSz^^^ 


■2w 


B^<p  ay 

dzdt  dt^ 


=  0 


(8) 


where 


u 


B<p 

dx' 


(9) 


Once  Eq.  (8)  is  solved  for  (f>,  the  Bernoulli -Kelvin  equation 
(4),  along  with  Eqs.  (5)  and  (6),  can  be  used  to  determine  the 
pressure  F  as  a  function  of  the  velocity  potential  <f>.  Therefore, 
in  terms  of  <p,  the  pressure  coefficient  C,,  defined  by 


C,= 


yP^M^ 


(10) 


can  be  expressed  as 
2 


yMi 


where  Af„,  the  reference  Mach  number,  is  given  by 


(M) 


(12) 


Tradirioml  llodogripl)  Methods  and  Exact  Sohrttons 

Early  investigators  like  Molenbrock'  and  Chaplygin^  were 
able  to  show  that  if  the  steady  approximation  of  Eq.  (8)  were 
mapped  onto  the  hodograph  plane  using  the  Legendre  trans¬ 
formation,  it  is  possible  to  determine  the  exact  fundamental 
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solutions  (to  this  nonlinear  equation).  Of  course,  this  was  a 
breakthrough  in  the  study  of  transonic  flow.  However,  noth¬ 
ing  seemed  to  have  been  said  about  the  applicability  of  this 
type  of  transformation  for  solving  the  unsteady  problem.  In 
an  eflbrt  to  determine  if  the  three-dimensional  form  of  Legen¬ 
dre’s  transformation  can  be  used  for  solving  Eq.  (8),  the 
following  results  were  obtained 
Consider  the  transformation 

<f)  =  xC  +  zri  IX  -  X  (13) 


is  the  Legendre 

transformation 

in  three 

50 
^  d.x' 

,3 

II 

II 

sx 

dx 

dq 

Equation  (13)  transforms  Eq.  (8)  into  its  counterpart  in  the 
hodograph  space  given  by 

(c^  -  f -  x^,]  +  2(>ilx,arr  -  xi,x„] 

+  (c^  “  ■■  Z{t)  “  —  XnXfd 

-  2rilXi;X„  -  X(,X(ni  ~  Xc{X„n  +  Zf-i  =  ® 

Obviously,  by  comparing  Eqs.  (8)  and  ( 1 5)  it  can  be  seen  that 
Eq.  (15)  does  not  look  easier  than  Eq.  (8).  This  type  of 
comparison,  therefore,  explains  why  the  traditional  approach 
is  not  effective  for  solving  Eq.  (8). 

New  Hodogriph  Mapptng  Scheme 

Before  formulating  the  new  mapping  scheme  (hat  can  trans¬ 
form  Eq.  (8)  into  its  relatively  easier  counterpart  in  the 
hodograph  plane,  a  nondimensionalization  scheme,  which 
uses  the  following  set  of  affine  transformations,  is  considered. 

Afflne  Transfonoatkim 

bio 

=  •*=<'*0.  x  =  bzo,  l  =  —  (16) 

Equation  (16)  transforms  Eq.  (8)  into  its  nondimensionalized 
counterpart  in  the  affine  space,  given  by 


d^(bn  ^^00  J  ^^00 

,  3'0O 


d^0o  b^<t>o  5^00  _  j.  / , 

“’’dxgdlo  "’’’bzo^io  ^<0 

Although  b  in  Eqs.  (16)  can  be  any  arbitrary  length  quantity, 
it  is  convenient,  at  least  for  computational  purposes,  to 
choose  b  as  the  chord  length  of  a  wing  section. 

Equation  (17)  is,  therefore,  the  nonlinear  velocity  potential 
equation  in  the  affine  space,  where 


500 

5.to' 

500 

’*’""5r„’ 

11 

(18a) 

II 

II 

e 

c 

Co  =  — 

9*- 

(18b) 

Hodogriph  Trinfonnillon 

Consider  to  be  a  harmonic  function  of  the  form  given  by 
it>o(Xg,Zo.io)  =  nn-iiklg^nn-t 

=  <!>„( t„,r„),  mt  +  t  ^  kfo  ^  (n  +  I  )n  -  «  (19) 


where  t  -»0,  n  =0,1, 2, 3,....  This  is  essentially  a  quasi-unsteady 
assumption.  Therefore,  it  can  be  argued  that  this  analysis  is 
not  really  unsteady.  Clearly  a  response  to  that  argument  is 
that  the  analysis  is  probably  the  first  one  in  almost  one 
century  to  show  that  transonic  (nonlinear)  flows  other  than 
two-dimensional  or  unsteady  one-dimensional  flows  can  be 
studied  in  the  hodograph  space.  It  thus  provides  an  aspect  of 
transonic  flow  that  has  never  been  seen  before.  If  this  assump¬ 
tion  seems  crazy  to  some  readers,  it  should  be  noted  that  it  is 
not  without  precedent.  For  example,  this  assumption  is  not 
likely  to  be  crazier  than  the  famous  von  Karman-Tsien 
assumptions  that  led  to  the  definition  of  a  gas  known  as 
“Katman-Tsien”  gas.  This  gas  has  the  ratio  of  the  specific 
heats  to  be  equal  to  —  I,  i.e.,  y  =  —  I.  However,  the  Karman- 
Tsien’s  approximation  did  provide  good  benchmark  results 
that  compared  favorably  even  with  experimental  results.  The 
analysis  received  the  blessing  of  such  important  scientists  as 
Sir  M.  J.  Lighthill — particularly  pp.  367-373  of  Ref.  28. 

A  second  example  is  the  fact  that  the  so-called  small-distur¬ 
bance  theory  approximation  has  been  known  to  give  good 
results.  Therefore,  it  would  seem  that  the  primary  (most 
important)  thing  in  any  transonic  flow  analysis  is  the  preser¬ 
vation  of  the  flow's  nonlinear  character.  Every  other  thing 
seems  to  be  secondary. 

Furthermore,  it  should  be  pointed  out  that  as  a  conse¬ 
quence  of  the  assumptions  in  Eq.  (19),  if  the  steady-state 
solution  0o(xo.zo)  shocks,  such  shocks  would  be  retained 
in  the  quasi-unsteady  solutions  as  well. 

Finally,  perhaps  it  also  should  be  remembered  that  all  of 
the  transonic  flow  equations  (and,  indeed,  all  mathematical 
physics  equations),  including  the  Navier  Stokes  equations,  are 
all  approximations. 

The  nondimensionalized  frequency  (or  Sirouhal  number)  k 
is,  as  a  result  of  Eqs.  (16),  given  by 


where  w  is  the  circular  frequency  of  oscillations.  If  Eqs.  (19) 
and  (20)  are  substituted  into  Eqs.  (17),  the  following  equation 
is  obtained; 


OXn  OZq  CZ  0 


=  lik(ul-<rwl)  (21 : 

Equation  (21)  also  can  be  written  in  the  following  manner: 


(cg-«o)^-2u„Ho-^ 


-t-(cS-w?)-^Vk^0„ 


=  2ikql 

where  qg,  the  resultant  velocity,  is  given  by 
9o  =  “o  +  **'0 

Now  consider  the  following  transformation: 

Xo  =  "o to  +  "’o^o  -  00 


bu' 


Xo  =  iff 


When  Eqs.  (24-26)  are  substituted  into  Eq.  (22),  the 
following  equation  in  the  hodograph  plane  must  be  solved  to 
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obtain  the  transformed  velocity  potential; 

(C^  -  «  )  JTj  +  2uw  ——  +  (c"  -  ivO  — j 
ow^  oudw  cu^ 


and  c  and  J  are  the  nondimensionalized  steady-state  velocity 
of  sound  and  Jacobian  of  the  transformation,  represented  by 
Eqs.  (24-26),  respectively. 

From  Eq.  (27),  it  is  seen  that  if  the  Jacobian  J  is  prescribed, 
ab  initio,  the  resulting  linear  equation  can  be  solved  closed- 
form  to  determine  the  exact  fundamental  solutions  for  the 
transformed  potential  x  ^nd,  hence,  0.  Consequently,  these 
solutions  can  be  pieced  together  to  determine  the  shape  of  an 
airfoil.  Notice  that  the  familiar  steady  hodograph  equation 
can  be  recovered  from  Eq.  (27)  if  k  is  set  equal  to  zero.  It  also 
should  be  noted  that  the  presence  of  the  Jacobian  in  Eq.  (27) 
should  not  necessarily  be  considered  as  an  added  problem  in 
comparison  to  the  steady  hodograph  equation  in  which  the 
Jacobian  does  not  appear  explicitly,  since  its  behavior  must  be 
studied  in  both  cases  in  order  to  ensure  unique  solutions. 

Solution  Methods 

Although  Morawetz^’  has  shown  that  continuous  solutions 
for  Eq.  (27)  (steady-state  approximation  was  used  to  arrive  at 
this  conclusion)  for  a  closed  body  in  the  transonic  regime  do 
not  exist,  investigators  like  Nieuwland’  have  shown  that 
fundamental  closed-form  solutions  of  the  steady  approxima¬ 
tion  of  Eq.  (27)  can  be  pieced  together  employing  the  appro¬ 
priate  boundary  conditions  to  obtain  interesting  transonic 
flows  over  airfoils.  In  this  paper  an  effort  is  made  to  obtain 
closed-form  fundamental  solutions  of  Eq.  27,  which  are  con¬ 
sequently  pieced  together  after  the  appropriate  boundary 
conditions  are  enforced  to  study  some  steady  and  unsteady 
transonic  flows. 

Polar  Coordinatca 

In  an  attempt  to  obtain  fundamental  solutions  of  Eq.  (27), 
it  is  helpful  to  transform  this  hodograph  equation  of  motion 
into  a  polar  coordinate  system. 

Consider  the  following  transformations: 


u  =q  cosO,  w  =  q  sinfl 

(29) 

Equations  (29)  transform  Eq.  (27)  into  the  following  polar 
coordinate  counterpart: 

(30) 

where  J  is  now  the  Jacobian  of  the  transformation  into  the 
polar  coordinates  whose  relationship  to  J  can  be  seen  as 
follows; 

J  5(X,.2o) 

S{u,w) 

(31a) 

d(xo^o) 

S{q.0) 

(31b) 

d{x,2)  _  d{x„,Za)  d(q,0) 
d{u,w)  S(q,0)  f^fu.H’) 

(32) 

(3), 

q 

To  study  the  fundamental  solution  of  Eq.  (30),  consider  a 
flow  in  which 

J  =  J{g)  (34) 

Therefore,  a  general  solution  of  the  form 

Z  =  C(9)  cos(m9 -I- f)  (35) 

can  be  assumed  for  Eq.  (30)  in  order  to  obtain  the  following 
equation  for  Q{q). 


Equation  (36)  can  be  rewritten  in  a  slightly  different  manner 
if  the  following  definition  is  employed: 


where  M  is  the  local  Mach  number. 

Equation  (36),  in  light  of  Eq.  (37),  becomes 


Equation  (38)  must  be  solved  to  obtain  the  transformed 
velocity  potential  x  snd,  hence,  the  physical  space  velocity 
potential  <l>.  A  similar  approach  can  be  used  to  show  that  the 
following  equation  must  be  solved  to  obtain  the  transformed 
stream  function  given  by 

^  -  0(9)  cos(mff  +0  (39) 


-j^(l  -  + 0=0 


In  general,  Eqs.  (38)  and  (40)  are  hypergeometric  and, 
hence,  can  be  satisfied  by  combinations  of  power  and  loga¬ 
rithmic  series.  Finding  solutions  to  the  flow  around  a  closed 
body,  therefore,  becomes  dependent  on  piecing  these  types  of 
series  correctly.  Employing  such  fundamental  solutions  to 
construct  the  overall  solution  is  obviously  preferable  to  using 
arbitrary  series,  since  these  fundamental  solutions  are  “solid” 
in  that  they  satisfy  the  equation  of  motion. 

In  order  to  transform  Eqs.  (38)  and  (40)  into  the  familiar 
hypergeomelric  form,  consider  a  reference  q*  given  by 


and  a  change  of  independent  variable  given  by 


y  + 1  VflV 
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Equation  (42)  transfonns  Eq.  (40)  to  the  following 

An  exact  fundamental  solution  of  Eq.  (43)  can  be  given  by 
Q  =  B„i"''^F„(x)  (44a) 

where 

f  Jf)  cos(mO  + 1)  (44b) 

where  B„  are  arbitrary  constants  to  be  determined  by  the 
boundary  conditions  of  a  particular  flow  problem,  and  F„(T) 
satisfy  the  following  hypergeometric  equations; 


\  V  -  I  /  J 


+  ^  (m  +  l)(m  +  k^Ja)F„  =  0 


An  F„  satisfying  Eq.  (45)  is  a  hypergeometric  function 
given  by 


F„(t-)  =  f(d,h:J;f)  =  l4--j-^f 


d(d  +  nd(6  +0-7  , 
I  •  2  •  J(d  +  1)  * 


a 

d +  F=m -fi(l +k\);  ah=-^(m4-l) 

(m  +k^Jf,).  J=m  +  I  (47) 

In  terms  of  the  velocity  potential,  the  fundamental  solution 
is 

Z  =  -  B„r  -'^  ( I  -  t)  -''j^f„(f)  +  -^^"j^-jsin(m0  +  0  (48) 

Equations  (44)  and  (48)  represent  some  exact  fundamental 
solutions  for  the  stream  function  and  velocity  potential,  re¬ 
spectively,  in  the  hodograph  plane.  In  the  examples  shown  in 
this  paper,  the  stream  function  solutions  were  used  in  con¬ 
structing  the  flow  solutions. 

Coapatcd  Exanpla 

The  examples  computed  in  this  paper  consist  of  the  pres¬ 
sure  variations  with  Mach  numbers  and  Strouhal  numbers 
(reduced  frequencies)  for  some  chosen  points  on  a  flat  plate  in 
a  transonic  flow  and  what  may  be  considered  as  a  preliminary 
attempt  to  compute  the  steady  and  unsteady  transonic  pres¬ 
sure  distributions  around  a  70-10-13  supercritical  wing  section 
designed  in  the  early  seventies  by  Bauer,  Garabedian,  and 
Korn.‘  The  procedure  basically  consists  of  trying  to  construct 
compressible  (transonic)  flows  using  incompressible  flows 
around  an  elliptical  cylinder  resulting  in  a  shape  different 
from  the  cylinder  with  the  help  of  the  exact  fundamental 
potential  flow  solutions  above.  This  procedure  has  been  thor¬ 
oughly  documented  by  Nieuwland’  for  the  steady  flow  prob¬ 


lem.  Hence,  no  attempt  has  been  made  in  this  paper  to 
redocument  the  basis  of  this  approach.  Therefore,  readers  are 
referred  to  Ref.  5.  Suffice  it  to  say  that  the  approach  involves 
a  lot  of  “bookkeeping”  and  patience.  Furthermore,  it  must  be 
pointed  out  that  this  procedure  can  (and  should)  be  better 
automated.  This  endeavor  is  currently  in  progress. 

The  flat  plate  results  shown  in  Figs.  1-6  have  basically  two 
main  features: 

1)  Dips  in  the  pressure  distributions  when  the  Mach  num¬ 
bers  are  varied  in  the  transonic  regime  seem  to  agree  with  the 
general  trends  in  the  behavior  of  transonic  characteristics  such 
as  the  lift  coefficient,  the  so-called  transonic  dip,  which  is 
basically  the  loss  of  aeroelastic  stability  at  transonic  speeds. 

2)  The  pressure  distribution  at  a  particular  Mach  number 
at  a  particular  chord  length  seems  to  have  an  optimum 
reduced  frequency  (Strouhal  number). 


Fig.  I  Pressure  coeffldent  «s  Mich  nomben  for  ■  flat  plilc. 


Fig.  2  Pressure  coeflicienl  rs  Mach  nnniliers  for  a  flat 


Hg.  .X  Pressure  coeflicienl  rs  Mach  lumibers  for  a  flat  plite. 
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Fig.  4  Prcssore  cocfliciciit  n  rtdaced  fra|iiciicy  for  a  flat  plate. 


Fig.  6  Pr«*ar«  coetndeat  n  rcdaccfl  tngacocy  for  a  flat  pUlc. 


Marble”  has  shown  similar  trends  for  a  quasi-one-dimen- 
sional  flow.  The  result  shown  in  Fig.  7  for  a  70-10-13  super¬ 
critical  airfoil  designed  back  in  the  early  seventies,*  which  is 
preliminary  in  nature,  seems  to  indicate  the  feasibility  of  using 
the  approach  outlin^  in  this  paper  in  designing  unstwdy 
aerodynamic  characteristics  for  an  airfoil.  The  pressure  distri¬ 
bution  for  the  supercritical  airfoil.  70-10-13  shown  in  Fig.  7, 
is  for  a  freestream  Mach  number  ot  0.7  and  a  lift  coefficient 


Fig.  7  Prcssore  dtstribadoo  on  a  7t>-IO-l3  sapercritical  airfofl. 

(Cf^)  —  0.98,  and  the  airfoil  has  a  thickness  ratio  of  0.127.  The 
angle  of  attack  for  the  steady  computation  was  0  deg.  For  the 
unsteady  computations  it  is  assumed  that  kr,  =  0.1,  and  that 
the  airfoil  is  oscillating  about  the  half-chord  with  a  pitch 
amplitude  of  1  deg.  The  main  significance  of  Fig.  7  is  that  the 
approach  outlined  in  this  paper  can  be  used  to  accurately 
compute  the  pressure  distributions  for  supercritical  airfoils, 
since  these  pressure  distributions  compare  very  well  with  the 
results  obtained  in  Ref.  6  for  the  steady  flow  and  also 
compare  favorably  with  results  of  Ref.  IS  for  the  unsteady 
flow.  Because  of  the  fact  that  the  analysis  in  this  paper  is 
preliminary  in  nature,  a  more  general  set  of  conclusions  would 
have  to  await  more  refined  computations  that  are  anticipated 
in  the  near  future.  Finally,  it  must  be  pointed  out  that  the 
results  presented  in  this  paper  are  basically  shockless  invisdd 
results.  Therefore,  the  plots  shown  for  Mach  numbers  greater 
than  unity  can  only  ^  approximations  to  the  weak  shock 
solutions. 

Concluding  Remarks 

This  paper  has  attempted  to  establish  the  existence  of  exact 
closed-form  fundamental  solutions  to  the  two-dimensional 
nonlinear  unsteady  potential  aerodynamic  equations.  Evi¬ 
dence  indicating  why  the  traditional  hodograph  approach  is 
not  effective  for  solving  the  nonlinear  unsteady  two-dimen¬ 
sional  flow  equations  is  presented.  Therefore,  suitable  map¬ 
ping  functions  are  employed  to  transform  the  nonlinear 
potential  flow  equations  into  the  hodograph  plane.  An  exam¬ 
ination  of  the  transformed  flow  equations  in  the  hodogranh 
plane  reveals  that  if  the  Jacobian  of  the  transformation  is 
prescribed  ab  initio,  the  exact  closed-form  fundamental  solu¬ 
tions  for  the  velocity  potential  and  stream  functions  can  be 
obtained.  It  is  seen  that  such  Chaplygin  solutions  can  be  used 
in  conjunction  with  the  incompressible  flow  around  elliptical 
cylinders  to  construct  transonic  flows  over  interesting  shapes 
with  the  help  of  a  methodology  developed  several  decades  ago 
by  investigators  such  as  LighthilF  and  Nieuwland’  for  steady 
flows.  Computations  of  the  pressure  distributions  for  certain 
points  on  the  flat  plates  seem  to  indicate  that  dips  exist  in  the 
pressure  distributions  as  the  freestream  Mach  numbers  are 
varied  in  the  transonic  regime.  The  results  also  seem  to  show 
the  existence  of  optimum  reduced  frequencies  for  the  pressure 
distributions.  A  steady  and  an  unsteady  pressure  distribution 
are  also  computed  for  a  70-10-13  supercritical  airfoil  designed 
back  in  the  early  seventies.*  Although  the  approach  needs 
more  efficient  automation,  the  results  computed  show  that  it  is 
feasible  to  use  this  approach  to  solve  the  inverse  problem  of 
designing  airfoils  with  desired  unsteady  aerodynamic  charac¬ 
teristics.  This,  therefore,  implies  the  possibility  of  designing 
airfoils  that  satisfy  both  the  performance  and  stability  criteria 
simultaneously. 
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FORMULATION  OF  3-D  riODOGRAPH  METHOD  AND  SEPARABLE 
SOLUTIONS  FOR  NONLINEAR  TRANSONIC  FLOWS  t 

* 

Gabriel  Oyibo 
Polytechnic  University 
Farmingdale ,  New  York 

Abstract 

Formulation  of  a  three  dimensional  hodograph  technique  for 
transforming  the  full  nonlinear  potential  transonic  flow  equa¬ 
tion  from  the  physical  space  into  an  equivalent  linear  counter¬ 
part  in  the  hodograph  plane  is  presented.  A  very  careful  exam¬ 
ination  of  the  governing  nonlinear  equations  in  the  physical 
space  reveals  that  a  mild  constraint  on  the  energy  equation 
(which  may  even  enhance  the  accuracy  of  this  nonviscous  foripula- 
tion)  would  permit  the  separation  of  the  nonlinear  flow  equations 
for  an  aircraft  wing  into  a  sectional  component  and  spanwise  com- 
pnjaent.  This  separation  of  variables  normally  believed  to  be  pos¬ 
sible  only  for  linear  equations  seems  to  have  been  possible  (for 
the  non-linear  equations)  because  of  some  inherent  mathematical 
symmetry  of  the  3-D  nonlinear  flow  equations.  The  consequential 
"three-dimensional"  sectional  equation  is  eventually  transformed 
into  the  hodograph  plane  where  it  becomes  linear.  A  further 
transformation  of  these  linear  hodograph  equations  into  the  char¬ 
acteristic  hodograph  plane  provides  the  opportunity  of  obtaining 
the  nonlinear  flow  field  for  a  particular  set  of  boundary  condi¬ 
tions  by  just  solving  a  set  of  first  order  characteristic  equa- 


•Associate  Rrofcs&or,  Aerospace  Engr.  Dept.,  Associate  Fellow  AIAA 


tions.  The  necessary  computations  can  easily  be  carried  out  on 
a  computer.  Some  preliminary  computations  show  good  agreement 
with  previously  computed  data.  This  verification  therefore  pro¬ 
vides  confidence  that  the  new  tool  can  perhaps  be  used  in  an  in¬ 
verse  manner  to  design  a  new  family  of  3-D  lifting  surfaces  with 
great  potential. 
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orized  to  reproduce  and  distribute  reprints  for  governmental  purposes 
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Nomenclature 


a,b,d,fQ 


B, 


B- 


c 

P 


/Y 


C,C 


C^(i=l,2, .  .  .) 
E,  t 


X 

F 

M 


x,y ,  z 


-s 

v,w 


u,  w 


<pf<p 


X 


A,  a 


f,h 


q,  0 

r, q 


coefficients  of  hodograph  equation 

complex  constant  and  its  conjugate,  respectively 

pressure  coefficient  and  fraction  of  wing  half 

span  (zero  is  wing  root)  respectively 

speed  of  sound  and  2-D  component  of  3-D  speed  of 

sound,  respectively 

coefficients  of  separated  pliysical  equations 
ellipse  parameter  and  time,  respectively 
fraction  of  chord  length  (zero  is  leading  edge) 
term  in  hodograph  equation 
Mach  number 

Cartesian  coordinates  (flow,  spanwise,  and  verti¬ 
cal  directions,  respectively) 
flow  velocity  components 
2-D  components  of  3-D  flow  velocities 
velocity  potential  and  stream  function,  respec¬ 
tively 

transformed  potential 

section  component  velocity  potential  and  trans¬ 
formed  potential,  respectively 

spanwise  component  velocity  potential  and  separa¬ 
tion  constant,  respectively 
functions  in  the  hodograph  plane 
flow  speed,  and  flow  angle,  respectively 
ratio  of  specific  heats  and  maximum  speed,  respec 
tively 

characteristic  coordinates 
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p,k 


T,  X 


S,  t 


(  ) 


00 


air  density  and  circulation  parameter, 
respectively 

regular  solutions  in  the  complex  characteristic 
hodograph  plane 

hodograph  variable  and  characteristic  roots, 

‘  respectively 

conjugate  complex  characteristics 
freestream  quantities 
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INTRODUCTION 


In  the  late  1800  s,  certain  Investigators'-’  discovered  that  the  hodograpb 
transformation  can  be  a  very  powerful  tool  for  studying  nonlinear  Inviscld  fluid  flow 
problems.  Subsequently,  later  Investigators^"  ‘‘‘  were  able  to  reveal  that  this  method 
could  be  used  to  design  supercritical  airfoils  with  superior  transonic  performance 
characteristics.  These>  developments  constituted  a  breaktiirougii  In  the  study  of  transonic 
flow  problems,  which  are  characterized  by  equations  which  exhibit  strong  nonllnearitles. 
The  success  of  this  method  was  due  to  the  following  two  characteristics  of  solutions 
obtained  from  the  hodograph  equations  as  compared  to  solutions  obtained  from  the 
equations  in  physical  space.  i''lrst,  the  hodograph  equations  are  linear  allowing  soli^tlons  to 
be  superimposed,  which  Is  not  possible  In  physical  space  due  to  the  nonlinearity  of  the 
u,<>\ (.•riling  e(inations.  yecoiid.  the  poleiitlal  for  ronniil.'iling  an  Inverse  problem  In 
hodograph  space  thereby  allowing  the  desired  ilnld  How  Held  to  be  specified  as  Input  and 
leading  to  the  rc«|nired  geoiiietry  as  the  computed  output. 


Despite  the  sncces.s  of  this  vvoik,  use  of  the  hodograph  technique  has  unnecessarily 
snITered  from  a  iiutlor  drawback:  that  is.  that  there  Is  a  general  belief  that  the  hodograph 
..c'Chn((|ue  can  only  he  applied  to  tiro-dimcn.sjon.'d  .«-.i.e;»dy  flow  proldcms.  This  Implies  that 
neither  unsteady  two-diniPiisioiial  flow  nor  stea<ly  three-diinenslonal  problems  can  be 
studied  using  a  hodograph  approaeli  and  ilnis  iiaa  linihed  i|ie  iiiierest  In,  and  application 


of.  I  his  powerful  analysis  KmI. 


15 

However,  recent  research  performed  by  Oyibo  has  shown  that  the 
hodograph  transformation  is  not  limited  to  two-dimensional  steady 
flows.  The  work  reported  in  this  paper  is  a  direct  cosequence  of 
the  discovery  that,  contrary  to  the  aforementioned 

general  belief,  the  hodograph  technique  can  be  employed  to  trans¬ 
form  the  nonlinear  transonic  flow  problems  with  three  independent 
variables  into  an  equivalent  linear  set  of  equations  in  the  hodo¬ 
graph  plane. 

The  derivation  of  these  transformed  equations^marks  the  first 
time  in  fluid  mechanics  that  the  hodograph  technique 

has  been  extended  beyond  steady  two-dimensional  flow  using  full  po- 
^ntial  equations.  It  is  important  to  note  that  a  straight-forward 
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application  of  the  usual  hodograph  (Legendre's)  transformation  for 
the  case  of  three  independent  variables  leads  to  an  equation  (cf. 

quation  (4)  below)  which  is  even  more  nonlinear  than  its  counter¬ 
part  in  physical  space.  While  Equation  4  is  consistent  with  the 
earlier  research  work  published  by  Guderley^^  which  was  indepen- 
dently  verified  by  Cole  and  Cook  '  using  a  different  approach, 
small  disturbance  (not  full  potential)  equations  were  used.  The 
primary  achievement  of  the  new  method  was  to  find  an  alternate  way 
of  defining  the  transformation  so  as  to  result  in  a  linear  system 
of  equations.  The  resulting  equations  have  tremendous  significance 
for  the  solution  of  both  steady  and  unsteady  fluid  flow  problems  due 
to  a  number  of  features  of  the  hodograph  method  in  general  and  the 
new  equations  in  particular. 

I  li  'i.  '-iiicT  the  lrnii.9l'ornie<^l  crjiiaiioii.s  me  iiiip;ir.  Mipy  alk>w  superposition  of  solutions. 

'I'liis  iiiofuis  thill  any  coiiipleie  solnlioii  ran  bo  ronsiriirtofl  from  a  combination  of 
fiiinhinioiilal  solutioiis.  This  will  enable  an  analytic  determination  of  the  Influence  of  the 
iiipnl  parameters  on  the  solutions  and  hence  will  lead  to  a  better  understanding  of  the 

flow  physics.  This  is  not  possible  with  the  nonlinear  equations  in^ 
'.n  physical  space  where,  except  for  specialized  cases,  only  purely 
numerical  solutions  are  possible.  These  numerical  solutions  suffer 
from  the  drawback  that  they  are  both  approximate  and  that  it  is  dif¬ 
ficult  to  understand  .the  effect  that  the  important  physical  param¬ 
eters  have  on  the  solution  to  the  problem,  since  these  effects  have 
to  be  deduced  by  studying  the  results  from  a  number  of  different 
cases . 

Second,  as  shown  below  in  Section)  •  •  l-l.  ihe  form  of  the  newly  derived 

hodograph  equations  for  the  both  the  unsteady  2-D  case  and  the  steady  3-D  case  are 
similar  to  the  steady  2-D  hodograph  equations  In  that  the  highest  order  (second) 
derivative  terms  are  identical  for  all  three  cases.  The  difference  between  the  cases  shows 
up  in  the  Inclusion  of  lower  order  terms  which  are  not  present  In  the  steady  two- 
dimensional  iiodograph  equation.  Since  the  computational  solution  techniques  and  the 
characteristics  of  partial  differential  equations  are  usually  determined  by  the  form  of  the 
highest  derivative  terms,  and  since  that  form  Is  Identical  to  the  steady  2-D  hodograph 
equations,  then  the  solution  techniques  previously  developed  for  steady  two-dimensional 
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How  slioulcl  be  cllreclly  applicable  to  the  3-D  steady  case.  This  Is  particularly  Important  In 
he  case  of  transonic  flow  due  to  the  non-uniqueness  of  the  solutions  as  proven  by 

20  fi 

Morawetz  .  For  the  case  of  2-D  steady  flow,  Garabedlan  and  Korn  have  worked  out  a 
systematic  and  efficient  procedure  to  enable  "closed”  bodies  to  be  developed  in  a 
straightforward  manner.  Since  the  new  equations  are  so  similar  to  the  ones  solved  In  |6j, 
this  work  Is  directly  applicable  to  the  current  problem  and  forms  the  bases  for  excluding 

I 

non-reallstic  solutions  from  the  admissible  set. 

Finally,  and  perhaps  the  most  Important  point.  Is  that  the  transformed  equations  allow 
inverse  solutions  to  be  obtained,  In  which  the  pre.ssiire  field  Is  specified  and  the  body 
geometry  Is  calculated  from  closed  form  solutions.  From  the  point  of  view  of  design,  this 
is  the  desired  situation  rather  than  the  usual  trial  and  error  method  of  picking  the 
geometry  and  then  examining  the  resulting  calculated  flow  field.  Thus  solutions  obtained 
from  the  new  hodograph  method  could  lead  to  the  design  of  new  families  of  wing  shapes 
which  should  prove  to  have  significantly  reduced  shock  strengths  and  hence  lower  drag 
when  shocks  appear  at  off  design  conditions.  The  fuel  savings  for  the  commercial  airline 

industry  made  possible  by  such  an  efficient  wing  design  could  ap- 
roach  hundreds  of  millions  of  dollars  annually  and  thus  the  pro¬ 
posed  technique  as  outlined  in  this  paper  could  lead  to  immediate 
and  substantial  industrial  value. 

In  reference  15,  the  two-dimensional  unsteady  nonlinear  transon¬ 
ic  flow  case  was  treated  using  the  new  hodograph  method  and  shown 
to  lead  to  results  in  good  agreement  with  previous  computational 
solutions.  The  ease  of  obtaining  solutions  to  the  new  equations 
was  demonstrated  by  using  a  microcomputer  to  obtain  these  solu¬ 
tions.  This  should  allow  more  room  in  big  computers  like  the  Cray 
for  design  optimation  rather  than  just  solving  for  the  flow  field. 

Since  that  time,  the  derivation  of  the  governing  hodograph  equa¬ 
tions  in  three  variables  has  been  extended  to  the  important  case 
of  steady  3-D  flow.  It  is  the  purpose  of  this  proposed  paper  to 
delineate  the  formulation  of  the  3-D  hodograph  equations  and  to 
carry  out  a  preliminary  study  of  their  solutions  for  the  case  of 
steady  flow  around  lifting  surfaces  of  finite  span. 
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Equations  {f7?i4:)  jand  (l7«.)  probably  are  among  the  very  few  (perhaps 
the  first  set)  of  separable  nonlinear  equations.  It  can  be  seen 
that  as  a  conse<^uence  of  equations  either  of  the  following  con 
straints  on  the  energy  (or  Bernoulli)  equation  become  necessary. 


2  2 
^08  ^00 

(— ^  =  0  in  three  dimensions 


(^oo'C„)  =  (u^,cJA(a,y) 


(  7h) 


For  the  constraints  in  Equation  (7f)  to  be  enforced,  must 
approach  -l.  This  is  either  an  interesting  coincidence  or 
something  is  being  said  about  the  realities  of  the  isentropic 
assumptions  in  the  potential  flow  formulations.  This  is  because 
=  -1  is  the  Karman-Tsien  gas  which  has  been  shown  to  be  so 
accurate  that  the  results  using  this  gas  agrees  very  well  with 
experimental  results.  This  agreement  which  is  really  excellent, 
can  be  seen  in  Figure  F,5b  in  reference  19.  Figure  F,5b  further 
shows  clearly  that  Karman-Tsien  gas  is  very  accurate  even  for 
truly  transonic  flows  (flows  with  supersonic  bubbles).  The 
alternative  constraint  in  Equation  (7h)  could  be  shown  to  work 
provided  that  t{p,  represent  the  reference  section  free  stream 
velocity  and  speed  of  sound  respectively. 
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Tlieiefore,  for  sldatJy  3-D  flowa,  tlie  liodograpli  equaiioa  to  be  solved  is: 

-f  “  /1X5  “  /2X5;  +  /3X  =  «  ^  ®) 

.  .r  =  0,  Equation  (  -S)  reduces  to  the  rainlllar  2-D  hodograph  equation. 

Equation  (:8)  Is  similar  to  the  2-D  hodograph  equation  which  has  been  studied  for 
almost  a  century.  Therefore  any  methods  of  solution  for  the  2-D  hodograph  equations  are 
also  applicable  to  tills  new  equation.  Due  to  the  mathematical  symmetry  mentioned 
previously,  the  solution  takes  on  the  separated  form  given  by  Equation  (/^)  and  Equation 
( '  8)  becomes  the  2-D  sectional  analog  of  the  .3-D  flow  problem.  From  Equation  the 

following  form  are  feasible  solutions  for  the  spanwise  component: 

00 

/l(a.y)  =  77  /l„c  (  9) 

n  =0 

where  a„  are  constants  determined  by  the  boundary  conditions  of  the  flow  and  y  is  the 
spanwise  coordinate. 

The  general  solutions  of  the  sectional  component  of  Equation  (  8)  are  In  terms  of 
hypergeometrIc  series.  Ttiese  solutions  be  used  along  with  the  method  of  complex 

characteristics  In  which  the  flow  Is  mapped  into  the  unit  circle  In  the  characteristic 
Jiodograpli  plane  in  order  to  obtain  solutions  for  Equation  (  8).  The  goal  Is  then  to  obtain 
the  body  stream  function  that  encloses  the  particular  wing  section  (  e.g.  the  root  section). 
Thereafter  the  spanwise  component  is  combined  with  this  solution  In  accordance  with 
Equation  to  provide  the  flow  field  over  the  entire  3-D  lifting  surface. 

1.2.,  Nonlinear  Tliree-Diitiensioiiai  Transonic  Flow  Investigation  and 

Construction 

The  fundamental  exact  solutions  to  the  new  hodograph  equations  are  hypergeometrIc 
scries  In  the  hodograph  space.  As  mentioned  above,  the  equation  system  In  the  iiodograpli 
plane  closely  resembles  the  two-dimensional  steady  hodograph  case  for  which 
computational  techniques  are  well  established.  It  is  expected  that  th^e  solution 
techniques  will  therefore  be  directly  transferable  to  the  higher  dimensional  cases.  The 
techniques  that  have  been  successfully  applied  In  previous  hodograph  analyses  fall 
generally  Into  two  categories,  either  solutions  obtained  by  superposition  of  a  series  of 
fundamental  solutions  as  In  the  manner  of  Nleuwland*  and  as  was  also  done  In  Reference 
or  by  converting  the  problem  Into  an  Initial  v.alue  problem  (or  a 
c  icterlstlc  Initial  value  problem)  In'  using  the  method  of  complex  characteristics 
developed  by  (Jarabedian.'® 
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While  either  of  the  solution  techniques  can  be  used,  the  method  of 
complex  characteristics  has  more  appeal  since  the  location  of  sonic 
1 Ines  and  limiting  lines  appear  more  naturally  in  characteristic 
coordinates  and  thus  this  method  is  expected  to  be  more  promising 

in  this  ongoing i study .  However,  in  the  present  study  the  series 
method  way  used  for  constructing  the  flows  satisfying  the  necessary  ' 
boundary  conditions  for  3-D  wings  and  to  study  the  fundamental  solu¬ 
tions  carefully  to  determine  how  to  piece  them  together  properly  as 
was  done  in  Reference  15  for  the  2-D  unsteady  case.  The  details  of 
this  procedure  are  given  in  [15] .  An  alternate  method,  which  fol¬ 
lows  the  2-D  work  of  Garabedian,  is  expected  to  be  the  primary  analy¬ 
sis  tool  the  next  stage  of  the  study,  and  consists  of:  (i)  reducing  the 
linear  hodograph  equation  into  its  canonical  form,  (ii)  utilizing  the  meth 
od  of  complex  characteristics  to  map  the  particular  section  flow  on¬ 
to  the  unit  circle  in  the  characteristic  hodograph  plane,  (iii)  ex¬ 
amining  the  solutions  to  ensure  that  the  body  stream  function  en¬ 
closes  the  particular  wing  section. 

For  that  purpose,  Equation  (>«)  reduced  to  Its  canonical  form,  would  lead  to  the 
f'^'lowlng  system  of  first  order  equations: 

F-  F- 

=0  ,  +  -f  — w,,  =  0  ()0) 

c  a 

-  X_u;^  =0  ,  —  X^iTj  =0  (|1) 

-  zw^  -  xlTj  =  0  (12) 

where  ^  and  i]  are  the  complex  characteristics  and: 


(13) 


with: 


Notice  that  the  canonical  system  of  equations  are  first  order  equations 

representing  Equation  (18)  or  &)■  One  way  to  solve  tliese  equations  Is  to  map  the  flow 
field  In  the  physical  space  Into  the  region  Inside  a  unit  circle  In  the  complex  characteristic 
hodograph  plane  and  then  use  tlie  incon\presslble  solution  to  provide  Initial  data  for 

I 

solving  the  characteristic  Initial  value  problem  resulting  from  the  mapping.  This  approach 
(thoroughly  documented  In  |5|)  Is  efficiently  automated  for  the  steady  2-D  case  In  a 
computer  code  which  we  have  obtained  from  Prof.  Garabedlan.  We  have  extensively 
exercised  this  code  and  are  thoroughly  familiar  with  It,  and  hence  liave  a  major  part  of  the 
necessary  code  for  the  3-0  problem  In  hand  due  to  the  similarity  between  the  2-0  and  3-0 
equations.  A  full  run  of  this  code  takes  only  30  to  ^5  seconds  on  a  Cray  X-MP  when  run 
with  non-vectorized  code.  It  Is  expected  that  vectorizatloii  of  the  code  which  Is  currently 
underway  will  enable  a  fully  optimized  (shock-free)  3-D  wing  solu¬ 
tion  to  be  obtained  in  about  10  minutes  which  appears  to  be  signi¬ 
ficantly  faster  than  the  time  current  computational  transonic  3-D 
codes  would  require  to  obtain  an  optimal  flow. 


Investigation  of  Fundanieiital  Solution 


i-X-i 

in  general  either  of  the  two  methods  outlined  above  can  be  exam¬ 
ined  and  used.  In  that  process,  the  fundamental  solutions  can  be 
studied  very  carefully,  particularly  to  ensure  the  proper  handling 
of  the  singularities  lihe  the  limiting  lines,  which  are  the  seros 
of  the  Jacobian  of  the  transformation  and  the  branch  lines,  w  er 
the  Jacobian  is  infinite. 

It  is  recognised  that  though  the  spanwlse  and  sectional  compon¬ 
ent  of  the  flow  are  ,in  a  sense,  considered  separately,  the  fu  da- 
.ntal  problem  of  shochless  flow  is  basically  to  find  smooth  tran 

sonic  solutions  for  Eguatlon  ,1  ),  describing  a  3-D 

about  a  lifting  surface  of  finite  span.  However.  It  is  also 


20 

recognized  as  shown  by  Morawetz  #  that  smooth  transonic  solutions 
do  not  even  exist  for  all  sectional  (2-D)  shapes  for  the  lifting 
'rface  (let  alone  for  a  3-D  lifting  surface) .  The  new  hodograph 
transformation  helps  us  to  overcome  this  difficulty  by  making  it 
possible  to  separate  the  3-D  flow  into  the  two  components  mentioned 
above  and  solve  them  as  an  inverse  problem  instead.  If  the  span- 
wise  component  is  given  by  Equation  (9),  then  the  sectional  com¬ 
ponent  is  to  be  determined  from  Equations  (|0-/4) .  In  this  study 
therefore,  we  computed  smooth  transonic  flow  by  piecing  together  the 
hypergeometric  series  solutions  to  equations  18  as  in  reference  15 
from  a  desired  flow  field  and  find  the  finite  span  lifting  surface 
which  generates  it.  This  procedure  was  used  to  find  the  necessary 

a  and  hence  the  spanwise  variation  of  the  flow  field.  In  order  to 

\ 

compute  the  sectional  component  of  this  smooth  transonic  flow  by 
the  method  of  complex  characteristics,  the  variables  in  Equations 

can  be  extended  into  the  complex  domain  where  a  characteris¬ 
tic  initial  value  problem  along  the  complex  characteristics  may  be 
solved.  While  this  procedure  can  be  done  in  a  purely  analytic  man¬ 
ner,  it  can  also  be  carried  out  conveniently  on  a  computer.  In  the 
ture  work,  we  propose  to  use  the  computational  approach  by  modify¬ 
ing  the  existing  code  of  Garabedian  to  include  the  additional  terms 
appearing  in  the  sectional  equation  (Equation  '8)).  When  this  is 
done  as  described  below,  the  "3-D"  section  will  be  determined  as 
well  as  a  which  provides  the  unknown  in  Equation  ('<9)  for  the  span- 
wise  component  of  the  flow. 

Considering  Equation  (|3),  for  real  variables,  for  supersonic  flow,  and  X_  are  real 
and  Equation  (J8)  Is  hyperbolic.  This  means  that  the  Initial  value  problem,  defined  as 
specifying  X  and  z  on  any  non-characterlstic  curve  as  well  as  the  characteristic  Initial  value 
problem  defined  as  one  in  which  x  or  z  Is  specified  on  one  characteristic  of  each  family  are 
well  posed.  However,  for  subsonic  flow,  X^  and  X_  become  complex  conjugates  and 
Equation  (iS)  Is  elliptical.  Here,  both  the  Initial  and  the  characteristic  Initial  value 
problem  are  no  longer  well  posed  In  the  real  domain,  and  hence  boundary  values  are 
generally  needed  for  formulating  a  correctly  set  problem. 

CREDIBILITY  OF  THE  NEW  HODOGRAPH  METHOD 

As  has  been  mentioned  earlier,  the  2-D  unsteady  results  have  been  shown.  In  |15| 
(which  has  been  reviewed  and  accepted  for  publication  in  the  ALAA  Journal  } 

to  be  accurate  when  compared  with 

previously  published  purely  computational  results. 


For  the  3-D  steady  case,  h'lgure  1  sliows  a  comparison  of  preliminary  results  of  tlie 
calculation  of  the  flow  from  the  new  method  wltli  results  obtained  by  Sprelter'®  using  a 
local  linearization  technique  for  the  case  of  A/  =  i  flow  around  a  rectangular  wing  with  a 
parabolic  arc  section  and  an  aspect  ratio  of  7.01.  The  results  are  shown  to  be  In  good 
agreement,  with  the  differences  attributable  to  the  approximate  nature  of  Sprelter's 
analysis. 


hlgur^,,2,show  comparisons  between  the  results  obtained  using  the  new  method  and 
those  computed  using  the  FL022  computer  code  of  Jameson  and  Caughey*’.  This  code 
uses  a  finite  difference  technique  to  solve  the  transonic  full  potential  ^equations  in  non¬ 
conservative  form.  For  the  casesshown  of  a  rectangular  wing  at  .  the  results  are 

seen  to  be  In  Excellent  agreement  due  to  the  more  accurate  model  contained  In  the 
computer  code  than  in  the  previous  case.  We  thus  feel  that  these  preliminary  solutions 
give  credibility  to  the  3-D  analysis  and  give  confidence  as  to  Its  potential. 


i,^r£:levance  of  proposed  technique 

Critics  of  the  use,  for  design  purposes,  of  the  full  potential  equations  (  upon  which  the 
hodograph  method  Is  based  )  often  say  that  numerical  solutions  of  the  either  the  three- 
dimensional  Navier-Stokes  equations  or  the  Buler  equations  would  provide  a  more 
accurate  prediction  of  the  flow  field,  particularly  at  lower  Reynolds  numbers  or  In  cases 
Involving  strong  shocks.  However,  accurate  numerical  predictions  by  themselves  do  not 
necessarily  provide  design  solutions,  nor  do  they  provide  guidance  to  designers  ns  to  what 
to  do  about  poor  aerodynamic  performance  nor  any  other  adverse  conditions  resulting 
from  the  predicted  aerodynamic  characteristics,  in  contrast,  the  solutions  from  the 
hodograph  space  do  provide  these  capabilities.  In  addition.  It  Is  Important  to  keep  In  mind 
that  strong  shocks  are  highly  undesirable  from  a  design  point  of  view  and  should  be 
avoided.  Thus  the  capability  of  predicting  strong  shocks  In  the  flow  field  Is  not  Imperative 
for  design  point  calculations.  For  the  off-desIgn  condition.  It  Is  expected  that  for  the 
aerodynamic  shapes  developed  In  this  research,  much  weaker  shocks  are  expected  than  for 
shapes  designed  based  on  two  dimensional  analysis.  For  this  situation  solutions  to  the  full 
potential  equation  are  still  valid  and  the  simulation  of  shocks  developed  as  part  of  the 
proposed  work  should  be  adequate  to  predict  off  design  aerodynamic  performance. 

Note  that  It  Is  generally  believed  that  the  discovery  of  supercritical  wing  sections  was 
made  possible  by  the  use  of  the  hodograph  transformation  along  with  the  full  potential 
flow  equations  (although  classified  experimental  work  was  evolving  the  same  airfoil  shapes 
Independently).  As  a  result,  new  characteristic  families  of  aliToll  shapes  were  defined  and, 
j  conjunction  with  wind  tunnel  tests,  were  proven  to  be  the  types  of  shapes  that  can 
postpone  or  eliminate  local  shocks  at  transonic  speeds  (a  conceptual  revolution!). 
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In  retrospect,  It  was  not  necessary  to  use  a  more  accurate  theory  like  the  Navler- 
Stokes  equations  In  order  to  discover  such  an  Important.  Tundamental  and  revolutionary 

I 

esult.  In  fact,  It  has  been  shown  that  the  boundary  layer  correction  effects  can,  in  most 
cases,  be  satisfactorily  Incorporated  empirically  Into  the  Inviscid  formulation  (I.e.,  without 
necessarily  using  the  Navler-Stokes  equations).  In  addition,  the  difficulties  of  generating  an 
appropriate  grid  and  establishing  the  appropriate  far  field  boundary  conditions  and  their 
location  for  the  computational  solution  of  the  equations  still  are  problems  that  require 
considerable  numerical  experimentation. 

In  spite  of  the  successful  design  of  shock  free  airfoils  using  the  hodograph  method,  true 
shock  free  wings  did  not  evolve  from  this  work.  One  of  the  drawbacks  of  trying  to  use 
supercritical  airfoils  to  construct  3-D  lifting  surfaces  Is  the  Inherent  significant  sensitivity 
of  the  transonic  flow  field  to  changes  In  Its  physical  parameters  due  to  the  strong  non¬ 
linearity  of  the  equations.  For  example,  small  changes  In  a  given  set  of  transonic  now 
parameters  can  result  in  a  significantly  different  flow  field.  This  means,  for  Instance,  that 
a  supercritical  section  designed  to  be  shockless  at  a  Mach  number  of  0.85  may  end  up 
having  strong  shocks  at  a  Mach  number  of  0.8  or  0.9.  This  type  of  behavior  Is  very  typical 
of  a  nonlinear  system  (which  Is  wliat  the  transonic  flow  problem  Is).  The  principle  of 
superposition  which  Is  useful  In  linear  systems  does  not  apply  to  transonic  flow.  This, 
therefore,  explains  why  a  three-dimensional  wing,  designed  with  a  set  of  two-dimensional 
sections  along  Its  span  that  are  Individually  shockless,  does  not  necessarily  end  up  being  a 
shockless  wing.  Unfortunately  however,  because  the  hodograph  method  (the  principal  tool 
for  designing  shockless  sections)  has  been  thought  to  be  (up  until  now)  a  two-dimensional 
tool,  designers  have  had  no  choice  but  to  use  the  superposition  Idea  to  design  a  three- 
dimensional  wing;  this  linearization  practice  (superposition)  Is  fairly  accurate  In  most  flow 
regimes  except  the  transonic.  The  penalty  for  such  an  Inaccurate  practice  Is  the 
reappearance  of  strong  shocks  (and  consequent  deterioration  of  aerodynamic  performance). 
Therefore,  a  three-dimensional  hodograph  approach  should  not  only  eliminate  such  an 
incorrect  practice  but  should  also  provide  the  correct  method  for  distributing  the  sectional 
properties  (e.g..  thickness)  along  the  span  (providing  perhaps  new  characteristic  families  of 
wing  shapes)  that  may  provide  truly  shockless,  three-dimensional  transonic  wings.  It  Is 
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enipliaslzed  tliat  the  availability  of  a  design  procedure  based  on  the  three-dimensional 
liodograph  space  transformation  would  provide  Insight  not  available  from  any  other 
approach. 

The  treatment  of  shocks  Is  very  Important  for  transonic  flow.  When  the  hodograph 
approach  Is  used  as  an  Inverse  method,  the  shock  Is  not  treated,  since  the  goal  In  such  an 
approach  Is  to  produce  "shockless"  airfoils.  Some  Investigators  are  skeptical  about  the 
existence  of  "shockless"  airfoils,  mainly  because  of  the  possibility  of  the  appearance  of 
shocks  at  the  off-desIgn  points.  This  Is  a  genuine  concern  since  the  "shockless"  airfoils 
which  have  been  designed  using  a  two-dimensional  theory  have  to  operate  frequently  at 
off-desIgn  points  e.g.  operation  In  a  three-dimensional  or  unsteady  physical  environment 
causing  local  off  design  conditions  even  at  the  design  Mach  number.  The  proposed  work, 
which  deals  with  application  of  the  hodograph  technique  to  the  three-dimensional  wing, 
should  provide  design  points  that  are  more  practical  and  should  alleviate  such  concerns. 

Finally,  this  section  may  be  concluded  with  a  quotation  from  Sir  M.  J.  Llghthlll:^  "It 
seems  likely  that  any  general  theory  of  compressible  flow  applicable  to  problems  with 
regions  both  of  sub-  and  supersonic  flow  (such  problems  have  been  called  'trans-sonic') 
must  be  based  on  the  'hodograph  transformation'  (due  originally  to  Molenbroek  1890  and 
Chaplygin  1004)." 


1.5. . CONCLUSIONS 

in  this  paper  the  formulation  of  a  new  3-D  hodograph  technique 
is  presented.  It  is  shown  that  a  mild  constraint  on  the  energy/ 
Bernoulli  equation  permits  the  separation  of  the  nonlinear  3-D 
flow  equations  over  an  aircraft  wing  into  the  sectional  component 
and  spanwise  component.  These  two  component  equations  are  solved 
in  both  physical  space  (spanwise  component)  and  hodograph-plane 
(sectional  component)  and  used  for  constructing  flows  over  3-D 

wings . 


While  the  formulation  of  the  complex  characteristics  method  of 
solution  for  the  sectional  component  is  presented  in  this  paper, 
the  hypergeometric  series  solution  method  is  utilised  for  construct¬ 
ing  the  flows  whose  results  are  shown  in  this  preliminary  compute- 
.ion  wort.  A  good  agreement  is  seen  when  these  results  are  compared 
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(ktid.  I?/  «kJi. 

to  the  previous  work  of  Spreiter^  That  verification  of  the  pres¬ 
ent  work  provides  the  confidence  that  this  new  technique  which 
seems  to  have  significant  potential  also  seems  to  possess  some 
credibility. 

Acknowledgements :  The  author  acknowledges  the  computational  help 
provided  by  Prof.  J.  Bentson,  Mr.  John  Nutakor  and  Mr.  Robert 
Jesinov/ski  all  of  the  Aerospace  Engineering  Dept.  (Polytechnic 
University) .  He  also  Acknowledges  several  discussions  with  Profs 
Paul  Garabedian,  Julian  D.  Cole,  and  William  Sears  of  the  Courant 
Institute  of  Mathematical  Sciences  (NYU) ,  Rensselaer  Polytechnic 
Institute  and  University  of  Arizona,  respectively.  The  author  is 
also  grateful  to  Dr.  T.  C.  Tai  of  David  Taylor  Naval  Research  Cen¬ 
ter  (Washington)  for  his  interest  in  and  appreciation  of  the  new 
method  described  in  this  paper.  The  author  also  thanks  Mrs.  B. 
Hein  for  patiently  typing  the  manuscript. 

1.^  REFERENCES 

1.  P.  Moleiibroek.  "Uber  elnlge  Bewegungen  elnes  Gases  bel  Annahme  eincs 
Geschwln<JlgkeRspoteiitlales,".Arch.  Math.  Phys.,  Vol.  0,  pp.  157-105,  1800. 

2.  S.  A.  Chaplygin,  "On  Gas  Jets,"  Scl.  Mem.  Moscow  Un.  Math.  Phys.  No.  21,  pp.  1- 
21,  1001.  Transl.  NACA  TM  1003,  1014. 

3.  M.  J,  Lighthlll,  "The  Hodograph  Transformation  In  Transonic  Flow,  I.  Symmetric 
Channels.  II.  Auxiliary  Theorems  on  the  HypergeometrIc  Functions.  III.  Flow 
Round  a  Body,"  Proc.  R.  Soc.  London,  AlOl,  pp.  311-300,  1047. 

4.  K.  G.  Guderley,  "On  Transonic  Simplifications  of  the  Hodograph  Equation,"  Wright 
Air  Development  Center  Tech.,  Report  No.  53-183,  1053. 

5.  G.  Y.  Nleuwland,  "Transonic  Potential  Flow  Around  a  Family  of  Quasl-elllptical 
Aerofoil  Sections,"  N.L.R.  Report  TR  T172,  1007. 

0.  F.  Dauer,  P.  Garabedian,  and  D.  Korn,  "Supercritical  Wing  Sections,  Lecture  Notes 
In  Economics  and  Mathematical  Systems."  Vol.  00,  Springer- Verlag,  New  York,  1072. 

7,  J,  W.  Boerstoel,  "A  Survey  of  Symmetrical  Transonic  Potential  Flows  Around 
Quasl-elllptical  Aerofoil  Sections,"  NLR  TR  T130,  1007. 


8.  H.  Sobleczky,  "Related  Analytical,  Analog,  and  Numerical  Methods  In  Transonic 
AJrfoll  Design,"  AlAA  Paper  70-1550,  July  1079. 

0.  A.  Hassan,  H.  Sobleczky,  and  A.  H.  Seebass,  "Subsonic  Airfoils  with  a  Given 
Pressure  Distribution,"  ALAA  Journal,  Vol.  22,  No.  9,  pp.  1185-1101,  September 
1084.  ' 

•  10.  P.  Garabedlah,  "Partial  Differential  Equations,"- John  Wiley  &  Sons,  New  York, 

1984. 

11.  R.  Courant  and  D.  Hilbert,  "Methods  of  Mathematical  Physics," 
Vol.  II,  Interscience,  New  York,  1962. 

12.  H.  Bateman,  "Partial  Differential  Equations  of  Mathematical 
Physics,"  Dover  Publications,  New  York,  1944. 

13.  K.  G.  Guderley,  "The  Theory  of  Transonic  Flow,"  Pergamon 
Press,  Oxford,  1962. 

14.  I.  E.  Garrick  and  C.  Kaplan,  "On  the  Flow  of  a  Compressible 
Fluid  by  the  Hodograph  Method  II  -  Fundamental  Set  of  Particu¬ 
lar  Flow  Solutions  of  the  Chaplygin  Differential  Equation," 

NACA  Report  No.  790,  1944. 

15.  G.  A.  Oyibo,  "Exact  Closed-Form  Solutions  for  Nonlinear  Un¬ 
steady  Transonic  Aerodynamics,"  DGLR  Paper  No.  85-22.  Second 
Int'l  Symposium  on  Aeroelasticity  and  Structural  Dynamics, 
Technical  University  of  Aachen,  W.  Germany,  April  1-3,  1985. 
Accepted  for  publication  in  AlAA  J. 

16.  Guderley,  G.,  On  Transonic  Airfoil  Theory,  J.  Aero.  Sci.,  Oct. 
1956,  pp.  961-969. 

17.  Cole,  J.  D.,  Cook,  L.  P.,  Ziegler,  F.,  Finite  Span  Wings  at 
Sonic  Speed,  Mech.  Rech.  Comm.,  1(4),  1980,  pp.  253-260. 

18.  Spreiter,  J.  R. ,  "The  Local  Linearization  Method  in  Transonic 
Flow  Theory",  in  "Symposium  Trannssonicum"  Aachen,  ^ent.  2-T . I^L2 

19.  Sears,  W.  R. ,  "General  Theory  of  High-Speed  Aerodynamics," 
Princeton  University  Press,  Princeton,  1954. 

20.  C.  S.  Morawetz,  "On  the  Non-existence  of  Continuous  Transonic 
Flows  Past  Profiles,  I,"  Comm.  Pure  Appl.  Math.,  Vol.  9,  pp. 
45-68,  1956. 

21.  Cole,  J.  ’D.,  Cook,  L.  P.,  "Transonic  Aerodynamics"  North-Hol- 
land  Series  in  Applied  Mathematics  and  Mechanics,  1986. 


.•^FIGURES  ^ 


.  0.6  - 

i  -  NEW  HOOOCRAPH  RESULTS 

O  SPREITER’S  RESULTS  [IBl 


X/CM.O 


FIG.  I  PRESSURE  COEFFICIENT  VS  SPANWISE'  STATION 
M*I.O,  a>7.7  DEGREES.  RECTANGULAR  WING 
ASPECT  RATIO  *7.63.  67.  PARABOLIC  ARC  SECTION 


SUrLKCJUTICAL 

DIHENSIOIJAL 


WjNc;  DL^jU.M 


1  ili<LE 


HODOGRAPH  APPROACH 


Gabriel  a 
John  Nu'vzakor'  an'.l 
Pw.i  yk  '  cr  'I '  c 
Farni!'-:  .  :  Lc  , 


Kyuno-:.0(',  Jang' 


Significant 

responsible 


conr.ri;;'. 
tor  the 


cions  rrc: 
avolucioii 


1  var. 
of  tl 


is  ir.\' 
lio'.lco; 


applicable  ro  nonliivear  pocential  flov;.  Works  of 
and  ChaplygjrP  were  compi imenced  by  later  efforts 


A  ^  /6^ 

r/v  rci . 

Wvcigatcrs  have  been 
rc.r.h  nethodology 


Mo].enbroc]:* 
of  Linhthill’, 


Gudarloy^,  Cole  and  Oook'’^.  Recent  v  Ori;s 
Bauer,  Garabedian  an'l  Korn’,  Gai  aLee-.i 
Seebass’  established  the  applicability 
shockless  supercritical  secrion  desiun 


ijy  RorawLet a  ' ,  N iewland"  , 
ii^  n  ilai-.san  ,  Sobiecal'.y  and 
of  the  Tiiethodolcgy  to 
Credibility  of  this 


by 


methodology  has  since  bc?en  clearly  established  eyperimen.tally 
studies  such  as  W'hi tcomb‘' .  Jhis  narjred  a  ma  jor  breakthrough  in 
transonic  flow  gas  dynamics,  iiov.-ever,  duo  to  the  mathemctical  and 
perceptive  constraints,  this  powerful  netnodology  has  been  viewed 
as  restricted  to  2-D  flow  or  gas  dynamics.  Jicv/ever  Oyibo“'”,  in  a 
recent  finding  was  able  to  prove  that  j-D  hodograph  methods  exist 
and  can  be  used  to  study  2-D  unsteady  and  3-D  gas  dynamic  flows 
over  3-D  wings  of  finite  span.  It  was  further  shown  in 
references  12  and  13  that  it  is  possible  to  find  smooth  or 
uhockless  solutions  to  the  2-D  unsteady  and  3-D  gas  flow  around 
wings  sections  and  3-D  wings  of  finite  span  respectively.  The 
smooth,  shockfree  solutions  generated  in  12  and  13  were  generated 
using  hypergeometric  series  (winch  are  fundamental  solutions)  the 
singularities  of  which  are  forced  inside  the  wing  and  out  of  the 
flow  field.  Rererence  12  also  shov/ed  the  formulation  of  the  basic 
3-D  hodograph  equation  in  the  characteristic  hodograph  variables 
in  a  manner  similar  to  that  formulated  .by  Garabedian*  for  2-D 
wing  sections.  This  formulation  along  with  the  separability  of 
the  solution  (into  a  "sectional"  component  and  a  wing  spanv/ise 
component)  as  established  in  reference  makes  it  possible  to 
modify  and  use  the  2-D  methods  to  solve  the  sectional  problem 
combine  it  v;ith  the  spanv;ise  solution  to  construct  a  smooth  or 
shockfree  solution  around  rcali.stic  closed  3-D  finiue  wing 
shapes.  This  write  up,  therefore,  user  this  characteristic 
hodograph  formulation  from  reference  13 
the  2-D  code  from  the  courant  Institute 
for  designing  shockfree  3-D  wings. 

Figures  1~^  show  the  results 


as 

of 


and 


a  basis  for  modifying 
Mathematical  Sciences 


obtained  in  this  research.  The 
which  confirm  the  findincs  of  reference 
3-u  f:;n;v.e  recc.angular  wing  is 
wcwjr.  th.o  research  .seo^ns  to 
may  not  i.>c  the  only  way  to  delay 
.zc  shockv.'avo  effects. 


preliminary  conclusions  ( 

13)  seem  to  be  that  n  .shockfree 
possible  and  feasible.  In  other 
suggest  that  sweep  and/or  twist 
drag  rise  mach  numbers  or  minim 

*  Associate  Professor,  Aerospace  Engineering 
Associate  Fellow,  AIAA 

*  Graduate  Research  Fellows,  Aerospace  Fnqi nee.-ina  Department 


Department . 


124 


REFERENCES 


I.  . Molenbrook ,  P.,  "Uber  einige  Bewequngen  eines  Gases  bei 
Annahme  eines  Geschwindigkaitspotentiales , "  Arch.  Math.  Phys.. 
Yol.9,  1890,  pp.  157-195. 

2.  Chaplygin,  S.A.,  "On  Gas  Jets,"  Sci .  Mem  Moscow  Un. 

Hath. Phys.,  1904  Transl.  NACA  TM  1063,  No.  21,  1944,  pp.  1-21. 

3.  Lighthill,  M.J.,  "The  Hodograph  Transformation  in  Transonic 
Flow.  I.  Symmetric  Channels.  II  Auxiliary  Theorems  on  the 
Hypergeometric  Functions.  III.  Flow  around  a  Body,"  Proc.  R.  Soc. 
London,  A191,  1947,  pp  341-369. 

4.  Guderly,  K.G.,  "  On  Transonic  Simplification  of  the  Hodograph 
Equation,"  Wright  Air  Development  Center  Tech.,  Report  No.  53- 
183,  1953. 

5.  Morawetz,  C.S.,  "On  the  Non-Existence  of  Continuous  Transonic 
Flow  Past  Profiles,  I,"  Comm.  Pure  Appl .  Hath.,  Vol.  (,  1956,  pp 
45-68. 

6.  Niewland,  G.Y.,  "Transonic  Potential  Flow  Around  a  Family  of 
Quasi-elliptical  Aerofoil  Sections,"  N.L.R.  Report  TR  T172,  1967. 

7.  Garabedian,  P.,  "Partial  Differential  Equations,"  Wiley,  New 
.ork,  1964. 

8.  Bauer,  F. ,  Garabedian,  P. ,  and  Korn,  D. ,  "Supercritical  Wing 
Sections,  Lecture  Notes  in  Economics  and  Mathematical  Systems," 
Vol.  66,  Springer-Verlag ,  New  York,  1972. 

9.  Hassan,  h. ,  Sobiecczky,  H.  ,  and  Seebass ,  A.H.,  "Subsonic 
Airfoils  with  a  Given  Pressure  Distribution,"  AIAA  Journal, 

Vol. 22,  No.  9,  1984,  pp  1185-1191. 

10.  Cole,  J.D.,  and  Cook,  L.P.,  "Transonic  Aerodynamics,"  North- 
Holland  Series  in  Applied  Mathematics  and  Mechanics,  1986. 

II.  Whitcomb,  R.T. ,  "  A  Study  of  the  Zero-Lift  Drag-Rise 
Characteristics  of  Wing  Body  Combinations  near  Speed  of  Sound," 
NACA  TR  1273,  1956. 

12  Oyibo,  G.A.,  "  Exact  Closed-form  Solutions  for  Nonlinear 
Unsteady  Transonic  Aerodynamics,"  AIAA  Journal,  Vol.  27,  No.  li, 
1989,  pp  1572-157a. 

13.  Oyibo,  G.A. ,  "  Formulation  of  Three  Dimensional  Hodograph 
Method  and  Separable  Solutions  for  Nonlinear  Transonic  Flows," 
AIAA  Journal,  Vol.  28,  No.  10,  1990,  pp  1745-1750. 


FIGURE/  DISTRIBUTION  OF  MACH  NUMBER  ON  THE  UPPER  SURFACE  OF  A  WING 


FIGURE'Z  DISTRIBUTION  OF  NACH  NUMBER  ON  THE  LOWER  SURFACE  OF  A  WING 


IGURE 


IGURE  4-  LOWER  SURFACE 


0.J4  0  42  0  50  0.58  0.66 


FIGURES  ISO-MACH  ON  THE  UPPER  SURFACE  OF  WING 


0.27  0.35  0.43  0.51  0..58  0.66  0.74  0.82 


FIGURE^  ISO-MACH  ON  THE  LOWER  SURFAOE  OF  WING 


dgaiAjoN  HOviAi 


EXACT  CLOSED  FORM  SOLUTIONS  TO  THE  FULL  NAVIER  STOKES  EQUATIONS 
AND  NEW  PERCEPTIONS  FOR  FLUID  AND  GAS  DYNAMICS 

Gabriel  A-  Oyibo* 

Polytechnic  University 
Farmingdale,  New  York 

The  full  Navier  Stokes  Equations  are  carefully  investigated  using 
the  fundamental  concepts  of  transformation  group  theory.  This 
reveals  that  in  some  transformed  spaces  Navier  Stokes  equations 
can  become  very  tractable  mathematically  and  fairly  transparent 
physically.  The  careful  examination  and  exploitation  of  these 
findings  permitted  the  formulation  of  rather  general  set  of 
fundamental  exact  closed  form  solutions  to  these  equations 
generally  considered  to  be  the  most  sophisticated  mathematical 
description  for  continuum  fluid  and  gas  dynamics.  These 
fundamental  solutions  are  then  used  to  construct  the  solution  to 
the  2-D  steady  flow  over  a  cylinder,  a  problem  which  apparently 
has  never  enjoyed  exact  closed  form  solutions  using  the  full 
Navier  Stokes  Equations  in  the  history  of  fluid  or  gas  dynamics. 
For  incompressible  flow  with  constant  viscosity  the  laminar 
separation  points  are  predicted  and  are  found  to  be  in  agreement 
with  experimental  data.  Similarly  the  pressure  distribution  as 
well  as  the  drag  predicted  by  these  new  solutions  are  found  to  be 
in  very  good  agreement  with  experimental  data.  This  experimental 
data  verification  therefore  provides  the  necessary  basis  for 
confidence  in  these  new  solutions  and  their  legitimacy  to  be  used 
to  explore  and  evolve  new  perceptions  in  fluid  and  gas  dynamics. 

In  addition  they  should  serve  as  benchmark  solutions  for 
validating  the  CFD  codes. 
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Stokes  Solution  for  Pressure  Distribution 
of  a  Creeping  Flow  over  Cylinder 
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X  Experimental  results  (Schllchting) 

Fig.  4  - Present  results 
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Drag  Coefficient  for  a  cylinder  versus  Reynolds  Number 
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STATUS  OF  PUBLICATIONS  AND  GRADUATE  STUDIES 


The  significant  results  obtained  thus  far  in  this  research  program  are  being 
compiled  and  written  up  for  publication  in  suitable  technical  journals.  AIAA  Paper  No. 
86-1006,  a  paper  based  on  the  first  phase  of  research  has  been  prepared  and  was 
presented  at  the  AIAA/ASME/ASCE/AHS  27th  Structures,  Structural  Dynamics  and 
Materials  Conference  held  in  San  Antonio,  Texas  May  19-21,  1986.  A  paper,  "Exact 
Solutions  to  Aeroelastic  Oscillations  of  Composite  Aircraft  Wings  with  Warping 
Constraint  and  Elastic  Coupling"  has  been  accepted  for  publication  by  the  AIAA 
Journal.  Another  paper  entitled  "Some  Implications  of  the  Warping  Restraint  on  the 
Behavior  of  Composite  Anisotropic  Beams"  has  been  published  by  the  AIAA  Journal 
of  Aircraft.  In  addition,  a  paper  entitled  "Exact  Closed  Form  Solutions  for  Nonlinear 
Unsteady  Aerodynamics",  which  provides  the  basis  for  the  proposed  next  phase  of 
the  program,  has  been  published  by  the  AIAA  Journal  and  included  in  this  report.  In 
the  process  of  studying  unsteady  2-D  nonlinear  transonic  flows,  what  seems  to  be  a 
major  breakthrough  was  discovered:  this  breakthrough  is  the  discovery  that  the  #-D 
steady  nonlinear  transonic  flow  equations  can  be  transformed  into  linear  hodogragh 
equivalent.  This  discovery  essentially  reverses  a  hundred  years  state  of  thinking  and 
belief  in  the  scientific  and  mathematical  community  that  the  hodogragh  method  was 
limited  to  the  2-D  flows.  The  results  of  this  discovery  have  been  written  up  and 
published  in  the  AIAA  Journal.  In  addition  another  paper  AIAA  92-257  which  presents 
practical  shock  free  3-D  wing  for  transonic  flight  is  to  be  presented  in  June,  1 992  at 
the  AIAA  applied  Aerodynamics  Conference  in  Palo  Alto,  California.  This  research 
program  seems  to  be  seeing  the  beginning  of  what  seems  to  be  another  breakthrough. 


perhaps  even  more  significant  than  what  we  have  seen  thus  far.  This  is  what  seems 
to  be  a  discovery  that  transformations  methods  can  even  be  used  to  effectively  study 
the  ultimate  set  of  equations  in  continuum  gas  dynamics  known  as  the  Navier  Stokes 
set  of  equations.  This  new  finding  could  be  the  key  to  effectively  unlocking  the  secret 
of  gas  dynamics  and  aeroelastic  phenomena  for  real  gases  with  viscosity,  which  have 
been  blurry  at  best  thus  far.  This  finding  has  resulted  In  a  preliminary  paper  accepted 
to  be  presented  at  an  international  conference  to  be  held  in  Colorado  in  August  of 
1992.  It  is  being  proposed  that  the  investigation  of  these  new  findings  be  pursued  in 
a  proposed  next  phase  of  research.  This  is  because  if  the  findings  are  correct  this 
could  be  the  start  of  a  new  era  of  fluid  or  gas  dynamics  analysis. 

Mr.  John  Calleja  and  Mr.  George  Papadopoulos  have  finished  their  Masters  thesis' 
under  this  project.  The  support  of  AFOSR  under  Contract  F  449620-87-C-0046  and 
Grant89-0050  of  these  thesis'  are  aknowledged  in  the  thesis  writeups. 
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